Falsifying Field-based Dark Energy Models 

Genly Leon, Yoelsy Leyva, Emmanuel N. Saridakis, 
Osmel Martin and Rolando Cardenas 

December 2, 2009 



Department of Mathematics, 

Universidad Central de Las Villas, Santa Clara CP 54830, Cuba 
Department of Physics, 

Universidad Central de Las Villas, Santa Clara CP 54830, Cuba and 
Department of Physics, 

University of Athens, GR-15771 Athens, Greece 

Email: genly@uclv.edu. cu, yoelsy@uclv.edu. cu, msaridak@phys.uoa.gr, osmel@uclv.edu. cu 
and rcardenas@uclv.edu. cu 

Abstract 

We survey the application of specific tools to distinguish amongst the 
wide variety of dark energy models that are nowadays under investiga- 
tion. The first class of tools is more mathematical in character: the ap- 
plication of the theory of dynamical systems to select the better behaved 
models, with appropriate attractors in the past and future. The second 
class of tools is rather physical: the use of astrophysical observations to 
crack the degeneracy of classes of dark energy models. In this last case 
the observations related with structure formation are emphasized both 
in the linear and non-linear regimes. We exemplify several studies based 
on our research, such as quintom and quinstant dark energy ones. Quin- 
tom dark energy paradigm is a hybrid construction of quintessence and 
phantom fields, which does not suffer from fine-tuning problems associ- 
ated to phantom field and additionally it preserves the scaling behavior of 
quintessence. Quintom dark energy is motivated on theoretical grounds 
as an explanation for the crossing of the phantom divide, i.e. the smooth 
crossing of the dark energy state equation parameter below the value -1. 
On the other hand, quinstant dark energy is considered to be formed by 
quintessence and a negative cosmological constant, the inclusion of this 
later component allows for a viable mechanism to halt acceleration. We 
comment that the quinstant dark energy scenario gives good predictions 
for structure formation in the linear regime, but fails to do that in the 
non-linear one, for redshifts larger than one. We comment that there 
might still be some degree of arbitrariness in the selection of the best 
dark energy models. 
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1 Introduction 

The current accelerated expansion of our universe has been one of the most ac- 
tive fields in modern cosmology. Many cosmological models have been proposed 
to interpret this mysterious phenomenon, see e.g. [H [5] for recent reviews. The 
simplest candidate is a positive cosmological constant A [3j H] . It is well-known 
that its interpretation as the vacuum energy is problematic because of its exceed- 
ing smallness ^ . Notwithstanding its observational merits, the ACDM scenario 
is seriously plagued by the well known coincidence and fine tuning problems [B] 
which are the main motivations to look for alternative models. 

Dark energy (DE) models with two scalar fields (quintessence and phantom) 
have settled out explicitly and named quintom models [3 [H [9l [TOl [TTl [12l [131 
[IllIIlllIelEaEaiiaESlESlEaEHll^ The quintom paradigm is a hybrid 
construction of a quintessence component, usually modelled by a real scalar field 
that is minimally coupled to gravity, and a phantom field: a real scalar field 
-minimally coupled to gravity- with negative kinetic energy. Let us define the 
equation of state parameter of any cosmological fluid as it; = pressure/density. 
The simplest model of dark energy (vacuum energy or cosmological constant) is 
assumed to have w = —1. A key feature of quintom- like behavior is the cross- 
ing of the so called phantom divide, in which the equation of state parameter 
crosses through the value w = — l.LJ Quintom behavior (i.e., the w — —1 cross- 
ing) has been investigated in the context of h-essence cosmologies |14[ [TS] ; in 
the context of holographic dark energy [171 [H UHl EHl EI] ; inspired by string 
theory [22l |23l [24]; derived from spinor matter [25]; for arbitrary potentials 
[27l|28l|29l|30] ; using isomorphic models consisting of three coupled oscillators, 
one of which carries negative kinetic energy (particularly for investigating the 
dynamical behavior of massless quintom) |31j. The crossing of the phantom di- 
vide is also possible in the context of scalar tensor theories [32] [33l [34l [35l [36] 
as well as in modified theories of gravity . 

^In section|2]we refer briefly to observational evidence in favor the quintom DE model. 
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The cosmological evolution of quintom model with exponential potential 
has been examined, from the dynamical systems viewpoint, in |12| and |13[ 116] . 
The difference between [T2] and [T31 US] is that in the second case the potential 
considers the interaction between the conventional scalar field and the phantom 
field. In |T3] it had been proven that in the absence of interactions, the solution 
dominated by the phantom field should be the attractor of the system and the 
interaction does not affect its attractor behavior. In [16] the case in which 
the interaction term dominates against the mixed terms of the potential, was 
studied. It was proven there, that the hypothesis in [13 is correct only in the 
cases in which the existence of the phantom phase excludes the existence of 
scaling attractors (in which the energy density of the quintom field and the 
energy density of DM are proportional) . Some of this results were extended in 
[27], for arbitrary potentials. There it was settled down under what conditions 
on the potential it is possible to obtain scaling regimes. It was proved there, 
that for arbitrary potentials having asymptotic exponential behavior, scaling 
regimes are associated to the limit where the scalar fields diverge. Also it has 
been proven that the existence of phantom attractors in this framework is not 
generic and consequently the corresponding cosmological solutions lack the big 
rip singularity. 

In the first part of the chapter we investigate basic cosmological observables 
of quintom paradigm. We perform the cosmological perturbations analysis of 
quintom model for independent quadratic potentials. We investigate the evo- 
lution of quintom cosmology with exponential potentials in a background of a 
comoving perfect uid. First, we review the at FRW subcase (with dust back- 
ground). Then, we consider both negative and positive curvature FRW models. 
We construct two dynamical systems, one adapted to negative curvature and 
the other adapted to positive curvature. We characterize the critical points 
of the resulting systems. By devising well-defined monotonic functions we get 
global results for ever expanding and contracting models. We find the existence 
of orbits starting from and recoUapsing to a singularity (given by a massless 
scalar field cosmology) for positive curvature models. There is also a closed 
FRW solution with no scalar field starting from a big-bang and recollapsing 
to a "big-crunch" . We have determined conditions for the existence of differ- 
ent types of global attractors. Furthermore, our monotonic functions rule out 
periodic orbits, recurrent orbits or homoclinic orbits. We comment about the 
interplay between dynamical analysis and observational checking as tools for 
discriminate among different quintom proposals. 

A large variety of dark energy models suffers from the eternal acceleration 
problem, due to the exponential de-Sitter expansion. One of the consequences 
of the eternal acceleration is that a cosmic horizon appears (see e.g. [SS] for 
a further discussion). This problem is not strictly related to A, should we 
replace A with a quintessence scalar field, the universe should still be eternally 
accelerated finally reaching a de Sitter phase and hence again a finite cosmic 
horizon. In the second part of the chapter we explore from both observational 
testing and dynamical systems perspective a theoretical model to address the 
horizon problem. We consider an effective dark energy fluid as a source of the 
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accelerated expansion. We follow a model presented by some of us [351 131] whose 
dark energy component is the sum of a negative cosmological constant and a 
quintessence scalar field evolving under the action of an exponential potentiaj^. 
As a result, although the model is presently accelerating, eternal acceleration 
disappears and the universe ends in a Big Crunch like singularity in a finite 
time. Motivated by these theoretical virtues, we further explore this model from 
the observational point of view in order to see whether a negative A is indeed 
compatible with the astrophysical data at hand. We conclude that a negative 
A is indeed allowed and could represent a viable mechanism to halt eternal 
acceleration. We also explore the predictions of this class of model concerning 
the structure formation in the Universe. We conclude that this model give good 
predictions for structure formation in the linear regime, but fail to do so in the 
non-linear. 

2 Observational Evidence for Quintom Dark En- 
ergy Paradigm 

In this section we are going to refer briefly on the observational evidence that 
favor the quintom DE model. 

2.1 Basic observables 

In this subsection we examine the basic observational quantity, which is the 
dark energy (DE) Equation-of-State (EoS) parameter. In 2004, supernovae la 
data were accumulated, opening the road to constraint imposition on the time 
variation of DE EoS. In |44| uncorrelated and nearly model independent band 
power estimates (basing on the principal component analysis |45j ) of the EoS 
of DE and its density as a function of redshift were presented, by fitting to 
the SNIa data. Quite unexpectedly, they found marginal (2cr) evidence for 
■w{z) < — 1 at z < 0.2, which is consistent with other results in the literature 

gsmaiiiiiigiEoiiii]. 

The aforementioned result implied that the EoS of DE could indeed vary 
with time. Therefore, one could use a suitable parametrization of iude as a 
function of the redshift z, in order to satisfactory describe such a behavior. 
There are two well-studied paramctrizations. The first (ansatz A) is: 

wbe = wq + w'z , (1) 

where wq the DE EoS at present and w' an additional parameter. However, 
this parametrization is only valid at low redshift, since it suffers from severe 
divergences at high ones, for example ate the last scattering surface z ~ 1100. 
Therefore, a new, divergent-free ansatz (ansatz B) was proposed [5^ [55] : 

2; 

fi'DE ^ Wq + Wi{l - a) = Wo + Wi — — , (2) 

J- "]— z 

another point of view of the composite dark energy models can be found in | 40l 1411 1421 143| 

where the cosmon model is introduced 
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where a is the scale factor and wi = —dw/da. This parametrization exhibits a 
very good behavior at high redshifts. 

In [53] the authors used the "gold" sample of 157 SNIa, the low limit of 
cosmic ages and the HST prior, as well as the uniform weak prior on f2m/i^, to 
constrain the free parameters of above two DE parameterizations. As can be 




wq Wo 

Figure 1: (Color online) Two-dimensional contour-plots of the DE equation-of- 
state parameters, in two different parameterizations and using SNIa data. The 
left graph corresponds to ansatz A (expression (Qpj and the right graph are to 
ansatz B (expression From Ref. \54^ . 

seen in Fig[l]they found that the data seem to favor an evolving DE with the EoS 
being below —1 around the present epoch, while it was in the range w > — 1 in 
the near cosmological past. This result holds for both parametrizations ([l]),©, 
and in particular the best fit value of the EoS at present is wo < —1, while its 
"running" coefficient is larger than 0. 

Apart from the SNIa data, CMB and LSS data can be also used to study 
the variation of EoS of DE. In [3S], the authors used the first year WMAP, 
SDSS and 2dFGRS data to constrain different DE models. They indeed found 
that evidently the data favor a strongly time-dependent wde, and this result is 
consistent with similar project of the literature ^i5ll[5Hl[5ni[Sni[SIllSl[S3[Ml 
[55] . Using the latest 5-year WMAP data, combined with SNIa and BAO data, 
the constraints on the DE parameters of ansatz B are: wq = —1.06 ± 0.14 and 
wi = 0.36 ± 0.62 [68l [66l EZ] , and the corresponding contour plot is presented 
in FigEl 

In conclusion, as can be observed, the current observational data mildly 
favor WDE crossing the phantom divide during the evolution of universe. 
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Figure 2: (Color online) Two-dimensional contour-plot of the DE equation- 
of-state parameters, in parameterization ansatz B (expression (^), and using 
WMAP, BAO, SNIa data. From Ref. J6^. 



Let us make some comments here. First of all, we mention that the above 
results can also fit the basic ACDM paradigm, where dark energy is attributed 
to the simple cosmological constant. Thus, many authors believe that accord- 
ing to data resolution we can still trust the ACDM paradigm, and thus there is 
no need to introduced additional and more complex mechanisms. The second 
comment is the following: even if we accept that the results seem to favor a 
DE EoS below —1 at present, this does not necessarily means that a two field 
explanation (one canonical and one phantom, i.e the basic quintom model) is 
automatically justified. One can still result to < —1 through many differ- 
ent frameworks including modified gravity, braneworld constructions, stringy or 
strong-inspired models, spinor models etc [69'. Thus, in order to distinguish 
between these alternatives, one has to find more complicated signatures of the 
two-filed quintom model, apart from the simple observable of DE EoS. One step 
towards this direction is to investigate the perturbation spectrum of two-field 
quintom model, and then examine its relation to observations. 



7 



2.2 Perturbation theory and current observational con- 
straints 

In this subsection we study the perturbations of two-field quintom DE paradigm 
and the effects of these perturbations on the current observations. AdditionaUy, 
since it is important to check the consistency of this model at the classical level, 
it requires us to analyze the behavior of perturbations when the EoS crosses the 
cosmological constant boundary [70] , 



2.2.1 Analysis of perturbations in quintom cosmology 

In the following discussion on the quintom perturbations we will restrict our- 
selves to the two-field quintom model, with a Lagrangian: 

£ = £q + £p, (3) 

where 



CQ = -d^cj,,d^^<P^-V^{c^{) (4) 



1 

2 

describes the canonical (quintessence) component, and 



>Cp--i9^025''02- 1^2(02) (5) 

the phantom one. The equations of motion for the two scalar fields = 1,2) 
read 

dV 

^^ + 2H^^±a^—^ =0, (6) 

where the positive sign is for the quintessence and the minus sign for the phan- 
tom field. Although in general the two scalar fields could be coupled with each 
other, here for simplicity we neglect these interactions. 

Now, for a complete study on the perturbations, apart from the fluctuations 
of the fields, one has to consider also the metric perturbations. In the conformal 
Newtonian gauge the perturbed metric writes 

ds^ = a^(T)[(l + 2^)dT^ - (1 - 2^)dx'dx,]. (7) 

Using the notation of 71j, the perturbation equations satisfied by each of the 
quintom components are: 

s^ = -{\^wi){e,~?>^)-m(^ -w^8, , (8) 

k = -ni-3w.)9,--^9, + k'(^-f^S,~a, + ^] , (9) 

l+Wi \ 1+Wi J 

where 

= {ei^i)54>,, fi, = , (10) 



8 



and 



SP, = Sp, - 2V,'S^, ^Sp^ + ^ [3n{l - wf) + m] . (12) 



fc2 

Thus, combining Eqs. dS]), (O and we obtain 
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2ne, + — ^ — 5, + k'^ , (13) 

l + w. 



-(1 + w,){e, - 3$) - 37i:(i - w^)s, - m 



w, + m{i - wf) 



(14) 



Since the simple two-field quintom model is essentially a combination of a 
quintessence and a phantom field, one obtains the perturbation equations by 
combining the aforementioned equations. The corresponding variables for the 
quintom system are 

(15) 



J2i P^^^ 
P^ 

and 



^quintom 



(16) 



Note that for the quintessence component, — 1 < wi < 1, while for the phantom 
component, W2 < — 1- 

The two-field quintom model is characterized by the potentials Vi. Lets 
us not consider the simplified case of quadratic potentials Vi{<j>i) — ^mf(j)f. In 
general the perturbations of arise from two origins, namely from the adiabatic 
and the isocurvature modes. Using instead of Si the gauge invariant variable 
Q = — $ — and in addition the relation $ = in a universe without 

anisotropic stress, the equations HH) and can be rewritten as, 

H 

3 ^' V"' ' ^ ' 
e, = 2nd, + k\3C^ + '^'^) , (19) 



G = -i~c,[c^ + ^ + J:^o,] (18) 



where 

a = + 3H(1 - w,) = ao[ln(a'^|p, +k|)]. (20) 

In these expressions Q is the curvature perturbation on the uniform-density 
hypersurfaces for the i-component of the universe |72j . Usually, the isocurvature 
perturbations of (pi are characterized by the differences between the curvature 
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perturbation of the uniform-(/)i-density hypersurfaces and that of the uniform- 
radiation-density hypersurfaces, 



S.,r = 3(C^ - Cr), (21) 

where the subscript r stands for radiation. Here we assume there are no matter 
isocurvature perturbations, and thus — Cr- Ehminating Q in equations (jlSp 
and (jl9p . we obtain a second order equation for 9i, namely 

e, + (Q - 2n)ei + {dn ~2n + k^)6, = k^{A^ + c,$). (22) 

The general solutions of this inhomogencous differential equation, is the sum of 
the general solution of its homogeneous part with a special integration. In the 
following, we will show that the special integration corresponds to the adiabatic 
perturbation. 

As it assumed, before the era of DE domination, the universe was domi- 
nated by either radiation or dark matter. The perturbation equations for these 
background fluids read: 

C/ = -e//3 , 

^ -n{l - iwf)ef + k'^i^WfCf + {1 + 3u;/)$] . (23) 
From the Poisson equation 

a ^ ^ ^ 

on large scales we approximately acquire: 

Therefore, combining the equations above with H = 2/[(H-3w/)t], we get (note 
that numerically 9f ~ 0{k'^)(f) 

5 + 3wf 

^^ = -3(TT^'^ = ^°"^*" 

Of = '^^^±^^r . (26) 
^ 3(1 + ^ ^ 

Therefore, from (|22p we observe that there is a special solution which on large 
scales it is given approximately by 

Of = Of, (27) 

while (HH) leads to 

/•ad 



cr = c/- (28) 
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This indicates that the special integration of corresponds to the adiabatic 
perturbation. Hence, concerning the isocurvature perturbations of 0^, we can 
consider only the solution to the homogeneous part of 



9^ + {a - 2n)9, + (QH -2H + k^)di = . (29) 
These solutions are represented by 0,*'*° and Ct". The relation between them is 

cr = '-^-^ ■ (30) 

Since the general solution of Q is 

Q = + cr = Cr+ cr , m 

the isocurvature perturbations are simply Sir — 3Ci*°- 

In order to solve (|29|) . we need to know the forms of Ci and H as functions 
of time T. For this purpose, we solve the background equations ([6]). During the 
radiation dominated period, a = At , TL = 1/t and we thus we have 



and 



^-1/2 



^ 2\ , , , _2 



(32) 



(33) 



respectively, where A, Ai and Ai are constants, J,y{x) is the i^th order Bessel 
function and Ii,{x) is the vih order modified Bessel function. Since the masses 
are usually small in comparison with the expansion rate of the early universe 
rui <C 'H/a, we can approximate the (modified) Bessel functions as Ju{x) ~ 
x'^(ci +C2a;^) and Ii,{x) ^ x^{c.i +C2x'^). We mention that J_i/4 and /_i/4 are 
divergent when x ^ Q. Given these arguments we can see that it requires large 
initial values of (f)i and 4>i if A2 and A2 are not vanished. Imposing small initial 
values, which is the natural choice if the DE fields are assumed to survive after 
inflation, only Ai and Ai modes exist, so 0i will be proportional to in the 
leading order. Thus, the parameters Ci in ([^D]) will be Ci = lO/r (we have used 
\Pi -\-Pi\ — 0|/a^)- So, we acquire the solution of ((29|) . 



gi^so ^ ^-4[^^^ cog(fc^) ^ sin(fcr)]. (34) 

Therefore, 91'^° presents an oscillatory behavior, with an amplitude damping 
with the expansion of the universe. This fact leads the isocurvature perturba- 
tions C,Y° to decrease rapidly. If we choose large initial values for 0^ and i^j, A2 
and A2 modes are present, (t)i will be proportional to in the leading order 
and Ci — 0. Now the solution of (l29t is: 



= T[Dii cos(fcT) + Di2 sin(fcr)] . (35) 
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That is, 9f° will oscillate with an increasing amplitude, so CI**" remains constant 
on large scales. 

Similarly, during matter dominated era, a — Bt^ , H = 2/t, and thus the 
solutions for the fields 0.; respectively read 



Therefore, we do get the same conclusions with the analysis for the radiation 
dominated era. Firstly, choosing small initial values at the beginning of matter 
domination, we deduce that the isocurvature perturbations in (jji will decrease 
with time. On the contrary, for large initial values, the isocurvature perturba- 
tions remain constant on large scales. This behavior was expected, since in the 
case of large initial velocity the energy density of the scalar field is dominated by 
the kinetic term and it behaves like a fluid with w = 1 , and thus its isocurvature 
perturbation remains constant on large scales. However, on the other hand, the 
energy density of the scalar field will be dominated by the potential energy due 
to the slow rolling, that is it will behave like a cosmological constant and thus 
there are only tiny isocurvature perturbations in it. 

In summary, we have seen that the isocurvature perturbations in quintessence- 
like or phantom-like field under quadratical potentials decrease or remain con- 
stant at large scales, depending on the initial velocities. In other words, the 
isocurvature perturbations are stable on large scales, with their amplitude be- 
ing proportional to the value of Hubble parameter evaluated during the period 
of inflation Hinf (if indeed their quantum nature originates from inflation). In 
the case of a large Hinf, the isocurvature dark energy perturbations can be 
non-negligible and thus they will contribute to the observed CMB anisotropy 
[73l[74]. However, in the cases analyzed in this subsection, these isocurvature 
perturbations are negligible. Firstly, as mentioned above, large initial velocities 
are not possible if we desire the quintom fields to survive after inflation. Fur- 
thermore, even if the initial velocities are large at the beginning of the radiation 
domination, they will be reduced to a small value due to the small masses and 
the damping effect of Hubble expansion. 

In conclusion, we deduce that the contributions of DE isocurvature pertur- 
bations are not very large [75 J and thus for simplicity we assume that Hinf is 
small enough in order to make the isocurvature contributions negligible. There- 
fore, it is safe to focus only in the effects of the adiabatic perturbations of the 
quintom model. 

2.2.2 Signatures of perturbations in quintom scenario 

Let us now investigate the observational signatures of perturbations in quintom 
scenario. For this shake we use the perturbation equations (fl6|) and (flT)) . and we 




(36) 



and 




(37) 
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are based on the code of C AMB [75] . For simplicity we impose a flat geometry as 
a background, although this is not necessary. Moreover, we assume the fiducial 
background parameters to be — 0.042, f^DM = 0.231, ri_D_E = 0.727, where b 
stands for baryons, DM for dark matter and DE for dark energy, while today's 
Hubble constant is fixed at Hq = 69.255 km/s Mpc~^. We will calculate the 
efii'ects of perturbed quintom on CMB and LSS. 

In the two-field quintom model there are two parameters, namely the quintessence 
and phantom masses. When the quintessence mass is larger than the Hubble 
parameter, the field starts to oscillate and consequently one obtains an oscil- 
lating quintom. In the numerical analysis we will fix the phantom mass to be 
mp ~ 2.0 X 10~^'^Mpi, and we vary the quintessence mass with the typical values 



being mg 10^^°Mpi and 4 x 10^''°Mpi respectively. 



Oscillatory Quintom 

In Fig. [3]we depict the equation-of-state parameters as a function of the scale 
factor, for the aforementioned two parameter-sets, and additionally their cor- 
responding effects on observations. We clearly observe the quintom oscillating 
behavior as the mass of quintessence component increases. After reaching the 
w = —1 pivot for several times, w crosses —1 consequently with the phantom- 
component domination in dark energy. As a result, the quintom fields modifies 
the metric perturbations: (5goo = '^o.^^jSgu = 2a^$(5y and consequently they 
contribute to the late-time Integrated Sachs- Wolfe (ISW) effect. The ISW effect 
is an integrant of <i> 4- 5" over conformal time and wavenumber k. The above 
two specific quintom models yield quite different evolving $ -I- ^' as shown in 
the right panel of Fig. [3l where the scale is fc ~ 10"'^ Mpc~^. As we can see, 
the late time ISW effects differ significantly when DE perturbations are taken 
into account (solid lines) or not (dashed lines). 

ISW effects constitute an important part on large angular scales of CMB and 
on the matter power spectrum of LSS. For a constant EoS of phantom it has 
been shown that the low multipoles of CMB will get significantly enhanced when 
DE perturbations are neglected |77j . On the other hand for a matter-like scalar 
field, where the EoS is around zero, perturbations will also play an important 
role on the large scales of CMB [78]. Our results on CMB and LSS refiect the 
two combined effects of phantom and oscillating quintessence. We mention that 
while in the early studies of quintessence effects on CMB, one could consider a 
constant Weff instead: 

_ Jdafl{a)w{a) 
"^^^ = Jdania) ' 

this is not enough for the study of effects on SNIa, nor for CMB, when the EoS 
of DE has a very large variation with redshift, such as the model of oscillating 
quintom considered above. 

To analyze the oscillating quintom-model under the current observations, 
we perform a preliminary fitting to the first year WMAP TT and the TE 
temperature-polarization cross-power spectrum as well as the recently released 
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Figure 3: (Color online) Effects of the two-field oscillating quintom on the ob- 
servables. The phantom mass is fixed at 2.0 x 10~^'^Mpi and the quintessence 
mass at 10~^^ Mpi (thicker line) and 4.0 x 10~^°Mp/ (thinner line) respectively. 
The upper right graph depicts the evolution of the metric perturbations $ + 5* o/ 
the two models, with (solid lines) and without (dashed lines) DE perturbations. 
The scale is k ^ Mpc^^ . The lower left graph shows the CMB effects and 
the lower right one delineates the effects on the matter power- spectrum, with 
( solid lines ) and without ( dashed lines ) DE perturbations. From Ref. \70l . 
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157 "Gold" SNIa data [7S]. Following [SDl [HI] in all the fittings below we 
fix r = 0.17, ^mh'^ = 0.135 and flth^ = 0.022, setting the spectral index as 
715 — 0.95, and using the amplitude of the primordial spectrum as a continuous 
parameter. In the fittings of oscillating quintom we've fixed the phantom-mass 
to be mp ~ 6.2 x 10'^^ Mpi. Fig. ddehneates 3cr WMAP and SNIa constraints 
on the two-field quintom model, and in addition it shows the corresponding best 
fit values. The parameters mg and mp stand for the masses of quintessence and 
phantom respectively. In the left graph of FigUwe present the separate WMAP 
and SNIa constraints. The green(shaded) area is WMAP constraints on models 
where DE perturbations have been included, while the blue area (contour with 
solid lines) is the corresponding area without DE perturbations. The pertur- 
bations of DE have no effects on the geometric constraint of SNIa. The right 
graph shows the combined WMAP and SNIa constraints on the two-field quin- 
tom model with perturbations (green/shaded region) and without perturbations 
(red region/contour with solid lines). We conclude that the confidence regions 
indeed present a large difference, if the DE perturbations have been taken into 
account or not. 

Non-oscillatory Quintom 



As we have mentioned, the basic observables could also described by the 
simple cosmological constant. Thus, in order to distinguish the quintom model 
from the cosmological constant, we now consider a quintom scenario where w 
crosses —1 smoothly without oscillations. It is interesting to study the effects of 
this type of quintom model, with its effective EoS defined in ([38]) exactly equal 
to —1, on CMB and matter power spectrum. Indeed, we have realized such a 
quintom model in the lower right panel of Fig. O which can be easily given 
in the two-field model with a lighter quintessence mass. In this example we 
have set mg ~ 2.6 x 10~^^ Mpi,mp ~ 6.2 x 10~^^Mpj. Additionally, we assume 
that there is no initial kinetic energy. The initial value of the quintessence 
component is set to (jju = 0.226Mpi, while for the phantom part we impose 
(^2i — 6.64 X 10~^Mpi. We find that the EOS of quintom crosses — 1 at z ~ 0.15, 
which is consistent with the latest SNIa results. 

The model of quintom, which is mainly favored by current SNIa only, needs 
to be confronted with other observations in the framework of concordance cos- 
mology. Since SNIa offer the only direct detection of DE, this model is the 
most promising to be distinguished from the cosmological constant and other 
dynamical DE models which do not get across —1, by future SNIa projects on 
the low redshift (for illustrations see e.g. [llj). This is also the case for the 
quintom model in the full parameter space: it can be most directly tested in 
low redshift Type la supernova surveys. 

In the upper left panel of Fig. [5] we delineate the different ISW effects 
among the cosmological constant (red/light solid), the quintom model which 
gives Weff = — 1 with (blue/dark solid) and without (blue dashed) perturba- 
tions. Similarly to the previous oscillating case, the difference is very large 
when switching off quintom perturbations and much smaller when including the 
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Figure 4: (Color online) 3a WMAP and SNIa constraints on two-field quin- 
tom model, shown together with the best fit values, mq and mp denote the 
quintessence and phantom mass respectively. We have fixed mp ~ 6.2 X 
10~^^ Mpi and we have varied the value of mq. Left graph: separate WMAP 
and SNIa constraints. The green (shaded) area marks the WMAP constraints 
on models where DE perturbations have been included, while the blue area ( con- 
tour with solid lines ) corresponds to the case where DE perturbations have not 
been taken into account. Right graph: combined WMAP and SNIa constraints 
on the two-field quintom model with perturbations (green/shaded region) and 
without perturbations (red region/ contour with solid lines). From Ref. JTUjj. 
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Figure 5: (Color online) Comparison of the effects of the two-field quintom model 
withweff = —1 and of the simple cosmological constant, in CMB (WMAP), the 
metric perturbations $ + ^' (the scale is k ^ 10""^ Mpc^^) and the linear growth 
factor. The binned error bars in the upper right graph are WMAP TT and TE 
data JM From Ref {Wf . 
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perturbations. In the upper right panel we find that the quintom model can- 
not be distinguished from a cosmological constant in light of WMAP. The two 
models give almost exactly the same results in CMB TT and TE power spectra 
when including the perturbations. We deduce that the difference in CMB is 
hardly distinguishable even by cosmic variance. 



2.2.3 Breaking the degeneracy between quintom and cosmological 
constant scenarios 

So far we have see that CMB observations cannot distinguish between a quintom 
model with Weff = — 1 and a cosmological constant. Thus, in order to acquire 
distinctive signatures, we have to rely in other observations. To achieve that we 
need to consider the physical observables which can be affected by the evolving 
w sensitively. In comparison with the cosmological constant, such a quintom 
model exhibits a different evolution of the universe's expansion history, and in 
particular it gives rise to a different epoch of matter-radiation equality. The 
Hubble expansion parameter becomes: 

H=-^ = Ha[n„,a-^ + nra-^ + X]^/^ (39) 

where X, the energy density ratio of DE between the early times and today, is 
quite different between the quintom-CDM and ACDM. In the ACDM scenario 
X is simply a constant, while in general for DE models with varying energy 
density or EoS we obtain 

X ^QoEa-^e-^^'"'-''^'^^'"'. (40) 

Therefore, the two models will give different Hubble expansion rates. This is 
also the case between the quintom model with We// = — 1 in the left panel of 
Fig. [5] and a cosmological constant. 

Finally, we mention that different H leads directly to different behaviors of 
the growth factor. In particular, according to the linear perturbation theory all 
Fourier modes of the matter density perturbations grow at the same rate, that 
is the matter density perturbations are independent of k: 

Sk + nSk - 47rGa2pM4 = 0. (41) 

The growth factor Di{a) characterizes the growth of the matter density pertur- 
bations, namely Di{a) — 6k{a)/Sk{a = 1), and it is normalized to unity today. 
In the matter-dominated epoch we have Di{a) = a. Analytically Di{a) is often 
approximated by the Meszaros equation [Sj : 

5n^H{a) r da' 
'''^''^^^^L {a'H[a')/H,r ^^^^ 

Therefore, we can easily observe the difference between the quintom and cos- 
mological constant scenarios, due to the different Hubble expansion rates. In 
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particular, one needs to solve (|¥T|) numerically. In the lower left graph of Fig. 
[5] we show the difference of Di{a) between the quintom model with Weff — —1 
and the cosmological constant one. The difference in the linear growth func- 
tion is considerably large in the late time evolution and possibly distinguishable 
in future LSS surveys and in weak gravitational lensing (WGL) observations. 
WGL has emerged with a direct mapping of cosmic structures and it has been 
recently shown that the method of cosmic magnification tomography can be 
extremely efficient [53 |HS1 HZ] , which leaves a promising future for breaking the 
degeneracy between quintom and a cosmological constant. 



3 Exponential quintom: Phase space analysis 

In the following discussion on the quintom phase space analysis we restrict 
ourselves to the two-field quintom model, with a Lagrangian: 

^ = ^d^cj^d'-c^ - ^a^'/'SV - V{c^, (43) 

and we include, also, ordinary matter (a comoving perfect fluid) in the gravita- 
tional action. As in [16] we consider here the efective two-field potential 

V = "Voe"^'™'^+"'^\ (44) 

where the scalar field represents quintessence and ip represents a phantom 
field. For simplicity, we assume m > and n > 0. 

Quintom (non-conventional) cosmologies with exponential potentials has 
been investigated, from the dynamical systems approach, for instance, in refer- 
ences [TSdSKlB] (see section \3A\ for a brief review). 

We shall consider the Friedmann-Lcmaitre-Robertson- Walker (FLRW) line 
element: 

ds^ ^ -dt^ + a\t) (^^^^+a;2(d^2_^sin2^V)) , (45) 

where k = 1,0, —1, identifies the three types of FRW universes: closed, flat, and 
open, respectively. 

The field equations derived from are 



H = -H 



1 (^2 _ ^2^ _ 1^^^ _ 1^^ ^ (46) 

1(^2-^2^ -I- il/ff-ipM, (47) 

PM - -3HpM, (48) 

4> + 3H4> - VQmV = 0, (49) 

^ -I- -iHip + V6nV = 0, (50) 

where H = denotes de Hubble expansion scalar. 

The dot denotes derivative with respect the time t. We consider a pressureless 
perfect fluid (dust) as the background matter. 
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3.1 Flat FRW subcase 

To investigate the flat models we introduced the same normalized variables as 
in [IS]: {x^, x^, y), defined by 



_ _ 'f _ vv 

They are related through the Friedman equation [46] by x'^ — x'^ + y'^ = 1 — jffj < 
1. 

The dynamics in the space is given by the ordinary differential equations 

m- 



x'^ = l (3my2 + (q _ 2)x^) (52) 
x'„^~^{3ny^ -{q-2)x,^) (53) 



y' ^ ^{1 + q - 3{mx^ + nx^))y (54) 



defined in the phase space given by 

= {x= ix^,x^,y) ■.0<xl~xl+y^ < 1}. (55) 

Here the prime denotes differentiation with respect to a new time variable 
T = loga'^, where a is the scale factor of the space-time. The deceleration factor 
q = —da/a^ can then be written 

g = i (3 {xl - xl - y^) + 1) . (56) 

The system (|52j|54p admits the critical points points O, C±,T, P reported in 
[16]). In the table[l]the location, existence and deceleration factor of the critical 
points for m > 0, n > and y > 0- We use the notation 6 = — n^. 

By analyzing the sign of the real part of the normally-hyperbolic curves C± 
we get the following results (we are assuming m > and n > 0): 



• If m < n, C+ contains an infinite arc parameterized by x* such that 

* ^ —n—m\ 
''V ^ 5 



X* < " f^iig^i is a local source. C_ contains an infinite arc pa- 

rameterized by X* such that x* < that is a local source. 



• If m = 71, C-i- contains an infinite arc parameterized by a;* such that 
a;* < ^™ that is a local source. All of C_ is a local source. 

• If m > n there are two possibilities 

— If 5 < 1, all of C- is a local source. A finite arc of C+ parameterized 
by X* such that < 2:* < ~"+"^VT^ is a local source. 
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Table 1: Location, existence and deceleration factor of the critical points for 
m > 0, n > and y > 0. We use the notation 5 ~ m? — (from reference [16]). 



Name 






y 


Existence 


q 














All m and n 


1 
2 


c± 









All m and n 


2 


p 


m 


—n 


Vi-s 


5 < 1 


-1 + 36 


T 


m 
2S 


n 

~2S 


1 

2VS 


S>l/2 


1 
2 



— If (5 > 1, no part of C+ is a local source and all of C_ is a local source. 

Perhaps, the most appealing result in fl6j is that, by introducing properly 
defined monotonic functions and by making some numerical integrations, it was 
possible to identify heteroclinic sequences 

• Case i) For m < -^/n^ + 1/2, the point P is a stable node, whereas the 
point T does not exist. The heteroclinic sequence in this case is C± — > 
O — >P. 

• Case ii) For + 1/2 < m < \/rP- + 4/7, the point T is a stable node 
and the point P is a saddle. For these conditions the heteroclinic sequence 

is C± — yO — yT — > P. 

• Case iii) For ^Jn? + 4/7 < m < vT+"r?, the point T is a spiral node and 
the point P is a saddle. For these conditions the heteroclinic sequence is 
the same as in the former case. 

• Case iv) For m > vT+r? the point T is a spiral node whereas the point P 
does not exist. The heteroclinic sequence in this case is C_ — > O — > T. 

From the possibilities listed above, there still the possibility of other attrac- 
tors different from phantom ones in the exponential quintom scenario, particu- 
larly scaling attractors (T) . This fact is a counterexample of one of the result 
in [13]. However, we are aware on the small probability that T represents the 
actual stage on the universe evolution due this solution is matter-dominated. 
Another compelling result in quintom cosmology is that the existence of mono- 
tonic functions in the state space can rule out periodic orbits, homoclinic orbits, 
and other complex behaviour in invariant sets. If so, the dynamics is dominated 
by critical points (and possibly, heteroclinic orbits joining it). Additionally, 
some global results can be obtained. A similar approach (i.e., that of devising 
monotonic functions) for multiple scalar field cosmologies with matter was used 
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Table 2: Location and existence conditions for the critical points at infinity. 



Name 9i 62 Existence 

±f always 

P2 TT ±-| always 

f ±cos-i(-f) -7r<±cos-i(-f) <^,n^0 

P^ ^ ±cos^i(^) -TT < ±cos-i (f ) < 7r,n 7^ 

Pb 0\ < 61* < TT 

Pe 6*1 TT < 61* < TT 



in [Ml [891 [90] . However, in that work they do not consider phantom-like scalar 
fields, as we do here. 

3.1.1 Analysis at infinity. 

The numerical experiments in 16J suggest that there is an open set of orbits 
that tends to infinity. Let us investigate the dynamics at infinity. In order to 
do that we will use the central Poincare projection method. Thus, to obtain the 
critical points at infinity we introduce spherical coordinates (p is the inverse of 



r — + a;^ + y^, then, p ^ as r ^ (X)): 



X0 = i sin^i COS02, (57) 

y = i sin 01 sin 02, (58) 
P 

x^p — — cos 01 (59) 

where < 0i < tt and — tt < 02 < tt, and < p < 00. 

Defining the time derivative f = pf, the system (j52ll54p . can be written as 

p' = i (cos2 01 - cos(202) sin^ 0i) + 2n cos 0i sin^ 0i sin^ 02p + O (p^) . (60) 
and 



9[ = ncos(20i) sin 01 sin^ 02 - cos0i sin0i sin^ 02p + O {p^) , 

02 = (n cos 01 cos 02 + TO sin 01 ) sin 02 — cos 02 sin02p + O (p^) . (61) 
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Table 3: Stability of the critical points at infinity. We use the notation 5 = 
m? — -n? and = n cos 9* ±m sin 9* . 



Name 



(Ai,A2 



Stability 



P.' 



-^3 



P5 



{-n,n) 
{-n,n) 

V2S S 



V2S S_ 



(0,A+) 
(0,A-) 



> 

> 

> 0, 6 <0 

< 0, 6>0 

> 0, S<0 

< 0, 6>0 




< 



< 



Sir 



4 ^ ^1 ^ 4 

otherwise 

4 < ^1 < — 

otherwise 



saddle 
saddle 

source if n < 0, n < m < —n 
saddle otherwise 

source if n > 0, —n < m < n 
saddle otherwise 

nonhyperbolic 

nonhyperbolic 



Since equation (|60p does not depends of the radial component at the limit 
p — > 0, we can obtain the critical points at infinity by solving equations (j6ip in 
the limit p — > 0. Thus, the critical points at infinite must satisfy the compati- 
bility conditions 

cos(26'i)sine'i sin^ 6*2 = 0, 
(n cos 01 COS 02 + TO sin 01 ) sin 02 — 0. (62) 

First, we examine the stability of the pairs (0*, O2) satisfying the compatibil- 
ity conditions ([S^ in the plane 0i-02, and then, we examine the global stability 
by substituting in and analyzing the sign of p'(02, 02)- In table [5] it is offered 
information about the location and existence conditions of these critical points. 
In table [3] we summarize the stability properties of these critical points. 

Let us describe the cosmological solutions associated with the critical points 
at infinity. 

The cosmological solutions associated to the critical points P^ and P^ have 
the evolution rates (j)^ /V = 0, tp/ip — and H/ip = p/V6 ^ 0. H These so- 
lutions are always saddle points at infinity. The critical points P^ and P^ 
are sources provided n < 0, n < to < —n or n > 0, —n < m < n, respec- 
tively. They are saddle points otherwise. The associated cosmological solutions 

^Do not confuse p with the matter energy density, the latter denoted by pM- 
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to have the evolution rates 4''^/V — , = —m/n, and H/(f) = 

—np/{^/3m) 0, and H/ip = p/VS 0, whereas the associated cosmological 
solutions to have the evolution rates 4>'^/V = , (p/^p = —m/n, and 

iJ/^ EE np/(V3m) ^ 0, and H/^ ee -p/V3 ^ 0. The curves of critical points 
P5 and Pe are nonhyperbolic. The associated cosmological solutions have expan- 
sion rates (valid for 9^ ^ 7r/4) V/<j)'^ = 0,'j)/(p = tanO^, H/ip = psecOyVG ^ 0, 
and V/(p'^ = 0, (/)/(^ — — tanfl^, -ff/</5 = psec6'|/'\/6 ^ 0, respectively. 

3.2 Models with negative curvature 

In this section we investigate negative curvature models. 



3.2.1 Normalization, state space and dynamical system. 

For the investigation of negative curvature models we shall use the normalized 
variables: (x^, x,^, y, fl), defined by 



</> 'P VV PM 

This choice allows to recast the Friedmann equation as 

1 - {xl -xl+y^ + n)=nk>o, (64) 

where 

= fc = -l'0- (65) 

Thus, 

0<xl^xl+y^ + n<l. (66) 

Let us introduce the new time variable, r, such that r — s- —00 as i — s- and 
T — > +00 as t +00. Since the time direction must be preserved we can choose 
dr = 3eHdt where e = ±1 = sign(iJ). 

The field equations are 

= e (I {q-2) x^ + my^) , 
x'^ = e (I (<? - 2) a;^ - ny^) , 
y' = e (|(1 + q) - mx^, - nx^) y, 

n' ^ le{2q-l) n. (67) 

Where = 2 ^a;^ — — + ^fi, is the expression for the deceleration 
parameter. The DE EoS parameter, w, can be rewritten, in terms of the phase 
variables, as 

9 2 9 
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Table 4: Coordinates and existence conditions for the critical points of the 
systeni|67l We have used the notation S = m? — n? . The subindexes in the labels 
have the following meaning: the left subindex (denoted by e = ±1) indicates 
when the model is expanding (+)or contracting (— ); the right subindex denotes 
the sign of Xip (i.e., the sign of <j)) and it is displayed by the sign ±. 



Label 


Coordinates: 


Existence 




{xp,x^,y,n) 




±K± 


{±^l + xf,x*^,0,0) 


All m and n 


±M 


(0,0,0,0) 


All m and n 


±F 


(0,0,0,1) 


All m and n 



±SF {m,—n, eVT^, 0) 5 <1 

±MS (^,-^,i^,^rr^) 5>i 



Notice that the evolution equation (|67l c) is form invariant under the coor- 
dinate transformation y — > ey. Then, the sign of ey is invariant by proposition 
4.1 in [91], in such way that we can assume, without lost generality, for fixed e, 
ey > 0. Hence, for each choice of sign of e, the equations (p7)) define a flow in 
the phase space 

= {{xp, x^,y,n):0<xl^xl+y^ + n< 1, 

xl-xl + y^>0,fl>0,ey>0}. (69) 

3.2.2 Form invariance under coordinate trasformations. 

First recall that the positive "branch" (e = +1) describe the dynamics of models 
ever expanding and the negative "branch" (e = —1) describes the dynamics for 
contracting models. The system is form invariant under the change e — > — e, 
i.e., the system is symmetric under time- reversing. In this way it is enough to 
characterize de dynamics in 5'"'". 

3.2.3 Monotonia functions. 

Let be defined in the phase space (or ^f^, depending of the choice of e) the 
function 
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Table 5: DE EoS parameter (w), deceleration parameter (g), fractional energy 
densities, and eigenvalues of the perturbation matrix associated to the critical 

points of the system [ST] We use the notation A''^ = nx* ± rriyj 1 + x'l^'^.±K±, 
±F, ±SF ±MS corresponds to fc = 0, the eigenvalues of these points in the 
invariant set of zero-curvature models are the same as displayed in the table 
but the first from the left. 



Label w q rtm,^de,^k Eigenvalues 



± 



K± 1 2 0,1,0 |e,0,e(l- A±),e 



M - 0,0,1 -|e,-|e,ie ^■ 



±F - I 1,0,0 ie,-ie,-ie,ie 



±SF -1 + 25 -1 + 3(5 0,1,0 2 {d-l)e,{5-l)e,{S~l)e,{2S-l) 



±MS i 1-^7^,0 ^-e,-le,-Ue±J{-7+l] 



{nx,j, + mx<o)^ Vt^ , , , 

M = ^— ^ 3 , M' = -2eM. (70) 

(1 — a;0^ + x^'^ — — il) 

This is a monotonic function for > and nXfj^+mx^ ^ 0. Then, the existences 
of such monotonic function rule out periodic orbits, recurrent orbits, or homo- 
clinic orbits in the phase space and also, there is possible global results from the 
local stability analysis of critical points. Additionally, from the expresion of M 
one can see inmediatly that fl ^ 0, ot nxcf, + mxip ^ o \nx^ + mxip\ +oo 
(implying x^ or x^p or both diverge) or fi^ ^ asymptotically. 

3.2.4 Local analysis of critical points. 

By the discussion about the invariance of the system, it is sufficient characterize 
dynamically the critical points +K±, +M, +F +SF +CS y +MS, in the phase 
space 4'+. In tables [4] and [5l it is offered information about the location, exis- 
tence and eigenvalues of the critical points of the system ([57)) in the phase space 
((69|) (for each choice of e) and also, it is displayed the values of some cosmological 
parameters associated to the corresponding cosmological solutions. 

Now we shall investigate the local stability of the critical points (and curves 
of critical points). We shall characterize de associated cosmological solutions. 



26 



The set of critical points ±K± and the isolated critical points ±M are lo- 
cated in the invariant set of massless scalar field (MSF) cosmologies without 
matter. The isolated critical points ±F are located in the invariant set of MSF 
cosmologies with matter. 

The arcs of hyperbolae ±K± parameterized by the real value a;* denote 
cosmological models dominated by the energy density of DE {fide — * 1), partic- 
ularly by its kinetic energy. DE mimics a stiff fluid solution. Since this arc a set 
of critical points, then necessarily, they have a zero eigenvalue. They are local 
sources (and in general they constitute the past attractor in the phase space 



The isolated critical points ±M denote the Milne's universe. They are non- 
hyperbolic. The critical points ±F represent flat FRW solutions (dominated by 
matter). They are hyperbolic. For this points the quintom field vanishes, then, 
the DE's cosmological parameters are not applicable to this points. 

The stable manifold of _)_M is 3-dimensional and it is tangent at the point to 
the 3-dimensional space {x^,x^,fl) whereas the unstable one is 1-dimensional 
and tangent to the axis y. This means the the critical point _|_M is unstable to 
perturbations in y. The critical point have a 2-dimensional stable manifold 
tangent at the point to the plane {x^, x^) and a 2-dimensional unstable manifold 
tangent at the critical point to the plane (y, 17). 

The isolated critical points ±SF and ±CS denotes cosmological solutions 
dominated by quintom dark energy and curvature scaling solutions, respectively. 
These are located in the invariant set of MSF cosmologies without matter (O = 
0). The critical points ±MS (belonging to the invariant set of MSF cosmologies 
with matter {Vt > 0)) represent flat matter scaling solutions. 

The stable manifold of ^SF in ^f"*" is 4-dimensional provided 5 < 1/3. In this 
case j^SF is the global attractor on j^SF is a saddle with a 3-dimensional 
stable manifold, if | < (5 < ^ or 2-dimensional if | < 5 < 1. 

The isolated critical points ±CS are non- hyperbolic ii 5 = |. On the other 
hand, the critical points ±MS are non- hyperbolic if ^ = 5. 

+CS is stable (with a 4-dimensional stable manifold) and then, it is a global 
attractor provided | < <5 < | (in this case all the eigenvalues are real) or if 5 > | 
(in which case there exists two complex conjugated eigenvalues in such way that 
the orbits initially at the subspace spanned by the corresponding eigenvectors 
spiraling toward the critical point). 

Let us notice that +MS is the global attractor of the system (it have a 4- 
dimensional stable manifold) only ifO<7< |,^> 2 (where 7 denotes the 
barotropic index of the perfect fluid). Since we are assuming 7 = 1 (i.e., dust 
background) then, the critical point +MS is a saddle. It have a 3-dimensional 
stable manifold if ^ < ^ < | (in which case all the eigenvalues are real) or if 
(5 > I (in which case there are two complex conjugated eigenvalues and then the 
orbits initially at the subspace spanned by the corresponding eigenvalues spiral 
in towards the critical point). 



^+) provided na;* ± mJl + a;*^ < 1. 
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3.2.5 Bifurcations. 

Observe that the critical points ±MS and ±SF are the same as S — > . _|_ SF 
{_SF) coincide with a point in the arc {_K_) as S ^ 1~ . This values of 5 

where the critical points coincide correspond to bifurcations since the stability 
changes. 

3.2.6 Typical behavior. 

Once the attractors have been identified one can give a quantitative description 
of the physical behaviour of a typical open (fc = —1) quintom cosmology. For 
example, for ever expanding cosmologies, near the initial singularity the model 
behave as de flat FRW with stiff fluid (DE mimics a stiff fluid) represented 
by a critical point in or in depending on the selection of the free 

parameters m,n and a;^ (see table [B]). Whenever ^CS exists (i.e., provided 
(5 > ^) it is the global attractor of the system. In absence of this type of points, 
i.e., if (5 < i, the late time dynamics is determined by the critical point +SF, 
i.e., the universe will be accelerated, almost flat {flk — > 0) and dominated by DE 
(r^de !)• DE behaves like quintessence (—1 < g < 0, i.e., —1 < w < — i) or a 
phantom field {q < —1, i.e., w < —1) if (5 > or (5 < 0, respectively. This means 
that, typically, the ever expanding open quintom model crosses the phantom 
divide (DE EoS parameter have values less than —1). The intermediate 
dynamics will be governed by the critical points +C5, +MS, y +M, which have 
the highest lower-dimensional stable manifold. 

Table 6: Summary of attractors of the system [ST] Observe that, whenever 
exists, the solution dominated by curvature _CS {+CS) is the past (the future) 
attractor for e = — 1, i. e., for contracting models (e = 1, i.e., for expanding 

models). We use the notation A='= = nx* ± mJ 1 + x*'^. 



Restrictions 


Past attractor 


Future attractor 


e = -1 


-SF if^<5 
-CS if^>i 


-K±if A± > -1 


e = 1 


+K± ff A± < 1 


+5"^ if<5<i 
+CS if^>i 



For contracting models, the typical behavior, is in some way, the reverse of 

*For fiat models, is weii l^nown tlie, wlienever it exists, (i.e., provided 5 > i) tiie attractor 
is +AfS (denoted by T in |16)). Wlien we include curvature, the stability of the matter scaling 
solution is transferred to the curvature scaling solution, as we prove here. 
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the above. If (5 < | the early time dynamics is dominated by - CS. Otherwise, if 
5 > \, the past attractor is -SF, i.e., the model is accelerating, close to flatness 
(ilfc — > 0) and dominated by DE. The intermediate dynamics is dominated at 
large extent by the critical points -CS", ^MS, y -F, which have the highest 
lower-dimensional stable manifold. A typical model behaves at late times as a 
flat FRW universe with stiff fluid (i.e., ME mimics a stiff fluid) represented by 
the invariant sets or -K^, depending on the choice of the values of the 

free parameters m, n y . 

3.3 Models with positive curvature 

In this section we investigate positive curvature models we shall make use of the 
variables similar but not equal to those defined in [88] section VI. A. 

3.3.1 Normalization, state space and dynamical system. 

Let us introduce the normalization factor 



i) = 3x/iJ2 + a-2. (71) 

Observe that 

(i.e., at a singularity). This means that it is not possible that D vanishes at a 
finite time. 

Let us introduce the following normalized variables (Qo: x^, x^, y, Cl), given 

by 

Q„ = X0 ^^-^,x^ = ^-^,y=^,n^ ^. (72) 
From the Friedmann equation we find 



0<xl-xl + f = l-n<l (73) 

and by definition 

- 1 < Qo < 1- (74) 
By the restrictions ((73l [74]) . the state variables are in the state space 

* = {(Qo, X4„x^,y) : < xl - xl + f < 1,-1 < Qo < 1}. (75) 

As before, this state space is not compact. 
Let us introduce the time coordinate 

I _ d 3d 
~ df i) dt 

D has the evolution equation 

b' = -3Qob (xl -xl + ^n 
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where 

This equation decouples from the other evolution equations. Thus, a reduced 
set of evolution equations is obtained. 

x'^ ^3my'^ + 3Qox^ (-1 + S) , 

x'^ ^ -Sny"^ + 3Qox^ (-1 + S) , 

y' ^ (mx^ + nx^ - QqE) . (76) 

Where E = x^- xl + ^Cl. 

There is also an auxiliary evolution equation 

Cl' = -Qo [-2 {xl ~ xl) + (l - Cl)) Cl. (77) 

It is useful to express some cosmological parameters in terms of our state 
variables. El 

{n^,nde, nk,q) = {n,i-n,Ql-i,-i + 3s) /qI, 

and 

-9 -2 -2 

3.3.2 Invariance under coordinate transformations. 

Observe that the system ((76l [77l) is invariant under the transformation of coor- 
dinates 

{t, Qo,x^, x^,y, Cl^ (^-f , -Qo, -x^.ij, (ij . (78) 

Thus, it is sufficient to discuss the behaviour in one part of the phase space, 
the dynamics in the other part being obtained via the transformation (|78[) . In 
relation with the possible attractors of the system we will characterize those 
corresponding to the "positive" branch. The dynamical behavior of the critical 
points in the "negative" branch is determined by the transformation (|75|) . 



3.3.3 Monotonic functions. 

The function 

is monotonic in the regions Qq < and Qo > for Qq ^ 1, nxtf, + nx^p ^ 
0, r2 > 0. Hence, there can be no periodic orbits or recurrent orbits in the 

^We have defined Qi. = .h,.) = .,\..> . 
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interior of the phase space. Furthermore, it is possible to obtain global results. 
From the expression N we can immediately see that asymptotically Qq — > 1 or 
nx^ + mx^p ^ or — > 0. 

3.3.4 Local analysis of critical points. 

In the tables [7] and [H] it is summarized the location, existence conditions, some 
properties of the critical points and the eigenvalues of the linearized system 
around each critical point. 

Table 7: Critical points of the system (|76p . We use the same notation as in 
table H 



Label 


Coordinates: 

{Qo,X4„x^,y) 


Existence 


±K± 


{e,±^l + xf,x*p,0) 


All m and n 


±F 


(e,0,0,0) 


All m and n 


±SF 


(e, me, —ne, y/l — 6) 


d < 1 


±CS 




0<5<i 


±MS 


( m n \/r /i 1 \ 
^"'26' 25'2\/5'V 


5>i 



In the following we will characterize the dynamical behavior of the cosmo- 
logical solutions associated with them. 

The critical points ±K, ±F, ±SF and ±MS represents flat FRW solutions. 

The set of critical points +if± parameterised by the real value represents 
stiff fluid cosmological solutions (DE mimics a stiff fluid) . It is the past attrac- 

tor for ever expanding models provided nx* ± m^Jl + x^ < 1. As we proceed 

before, a simple application of the symmetry (|78|) . allows to the identification 
of the future attractor for collapsing models: the typical orbits tends asymptot- 
ically to -K± as f ^ cx) provided na;* ± 1 + x* ^ > —1, and — 1 < Qo < 0. 
This fact has interesting consequences. If x* is a fixed value and n and m are 
such that — 1 < na;* + 1 + ^ < 1, then, there exists one orbit of the type 

+K+ — > If n and m are such that —1 < nx* — 1 + ^ < 1, then, 

there is one orbit of the type +K- These are solutions starting from 

and recoUapsing to a singularity given by a MSF cosmology (see figure [6(b)| ) . 

The critical points ±F represent flat FRW solutions. They hyperbolic. For 
these points the scalar fields vanish, so the cosmological parameters associated 
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Table 8: DE EoS parameter (w), deceleration parameter (q), fractional energy 
densities, and eigenvalues of the perturbation matrix associated to the critical 
points of the system ([7^ . We use the notation A"*" na;* ± m^l + x^'^. When 

the flow is restricted to the invariant sets Qo = ±1, the eigenvalues associated 
to the critical points ±F, ±SF and ±MS and to the critical sets ±K±, are, in 
each case, the same as those displayed, but the first from the left. 



Label w q flmj^dc^k Eigenvalues 



± 



k± 1 2 0,1,0 4e,0,3(e- A±) ,3e 



±F - I 1,0,0 e,le,~le,~le 

±S> + -1 + 36 0,1,0 2 (35- l)e,3((5- l)e,3((5- l)e,3e 
±CS -i 0,^,1-^ -2V35e, - V35e ± V4 - 96, -VSSe 



±MS i 1-^,^,0 6±J(-7+| 



to DE are not applicable to these points. If (5 > |, the ustable (stable) manifold 
of +F {-F) is tangent to the critical point and parallel to the plane y — Qo- 
This means that there is an orbit connecting -^F and _F pointing towards -F 
in the direction of the Qo-axis. It represents the closed FRW solution with no 
scalar field starting from a big-bang at and recoUapsing to a "big-crunch" 
at (see figure 6(a)). 

The critical point ^SF represents a solution dominated by the scalar field 
(with non- vanishing potential). It can be the global attractor in the sets < 
Qo < 1 or Qo = 1 (i-e., for ever expanding models, or flat models) for the 
values of the parameters displayed in [9l It can be a phantom dominated solu- 
tion provided 6 < 0. It also can represent quintessence dominated or de Sitter 
solutions. 

The critical point exist if 6 > i. They represent flat matter scaling 

solutions, for which both the fluid and quintom are dynamically important. It 
is a saddle point. 

For < (5 < i there exists the critical points ±CS for which the matter 
is unimportant, but curvature is non- vanishing (Qq ^ 1) and tracks the scalar 
field. These are called curvature scaling solutions. The values of its cosmological 
parameters are the same as for ±CS (displayed in table [S]), but it represents a 
different cosmological solution with positive curvature. These critical points are 
typically saddle points. 

In table [HI where we present a summary of attractors for the quintom model 
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with k ^ 1. 



Table 9: Summary of attractors for the quintom model with k = 1 (system (1761) '). 
We use the notation A='= ~ nx* ± m^/ 1 + x*'^. 



Restrictions 


Past attractor Future attractor 




-SF if5<i 


-1< Qo < 


- 5 as above 


< Qo < 1 


+k± if A± < 1 +SF if ^ < 1 




+SF if(5<i 

as above ~ „ . ^ ^ , 
+MS if > 1 



3.3.5 Bifurcations. 

Observe that the critical points ±MS y ±5^^ coincides as 5 i^. ±CS y ±S~F 
coincides as (5 ^ ^ . Additionally, +SF {^SF) coincides with a point at the 
arc +K+ as (5 1^. For this values of 5 a bifurcation occurs. 

3.3.6 Typical behaviour. 

Once the attractors have been identified one can give a quantitative description 
of the physical behaviour of a typical closed quintom cosmology. For example, 
for ever expanding cosmologies, near the big-bang a typical model behaves like 
a flat FRW model with stiff fluid represented by the critical set +K+ or by 
+K-, depending on the choice of the values of the free parameters m,n and 
x*. If (5 < i and < Qo < 1 the late time dynamics is determined by +SF, 
(with the same physical properties as +SF). The intermediate dynamics will 
be governed to a large extent by the flxed points +CS, +MS, and +F, which 
have the highest lower-dimensional stable manifold. For flat models (i.e., in the 
invariant set Qo = 1), the late time dynamics is determined by the critical point 
+SF provided i5 < ^ or +MS provided S > ^. 

For contracting models, the typical behavior is, in one sense, the reverse of 
the above. If (5 < | and — 1 < Qo < the early time dynamics is determined by 
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(b) 

Figure 6: (Color online) T/ie collapse of quintom cosmologies with positive cur- 
vature for the values for the parameters m — 0.7, n = 0.3 and j = I. In (a) we 
have selected the initial conditions QoiO) — f^(0) — 0.9, and ^^(O) — x^{0) = 
y{0) — 0. The dashed (blue) line represents the evolution ofQo vs t (observe that 
Qa evolves from 1 to —1, and eventually takes zero value). The dotted (red) line 
represents vs t. This illustrates the existence of a closed FRW solution with 
no scalar field starting from a big-bang at -^F and recollapsing to a "big-crunch" 
at -F. In (b) we have selected the initial conditions Q{){Q) — 0.9, f2(0) = and 
x^{Q) = — ^/l + iip(O)^ with x^p{Q) = 0.3. The dashed (blue) line denotes Qq vs 
T. Observe that Qq goes from the value 1 to —1 (i.e., the model collapses). The 
dot- dashed line denotes the evolution of x'^ — x"^ vs t (which is identically equal 
to 1). The dotted (yellow) line denotes the value of Xtp vs r and the straight 
(red) line denotes the value of x^ vs r. This illustrates the existence of orbits of 
the type +K- — )■ By choosing X(f,{Q) = -^/l + ^^'(O)^, with the same initial 

conditions for the other variables, we obtain orbits of the type +Kj^ -^K^. 
These are solutions starting from and recollapsing to a singularity (given by a 
massless scalar field cosmology). 
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-SF. The intermediate dynamics will be governed to a large extent by the fixed 
points -CS, -MS, and -F, which have the highest lower-dimensional stable 
manifold. For flat models (i.e., in the invariant set Qq = —1), the early time 
dynamics is determined by the critical point -SF (or _MF) provided S < ^ 
(^ > A typical model behaves at late times like a flat FRW model with stifT 
fluid (i.e. the dark energy mimics a stiff fluid) represented by the critical set 
-K^ or by ^K^ depending on the choice of the values of the free parameters 
m, n and x* . 



4 Observational Evidence for Quinstant Dark 
Energy Paradigm 

4.1 The model 

Looking at the impressive amount of papers addressing the problem of cosmic 
acceleration clearly shows that two leading candidates to the dark energy throne 
are the old cosmological constant A and a scalar field evolving under the 
influence of its self - interaction potential V{(j)). 

In the usual approach, one adds either a scalar field or a cosmological con- 
stant term to the field equations. However, since what we see is only the final 
effect of the dark energy components, in principle nothing prevents us to add 
more than one single component provided that the effective dark energy fiuid 
coming out is able to explain the data at hand. Moreover, as we have hinted 
upon above, a single scalar field, while explaining cosmic speed up, leads to 
a problematic eternal acceleration. A possible way out of this problem has 
been proposed by some of us [351 ISSl [S2 through the introduction of a negative 
cosmological term. 

Motivated by those encouraging results, we therefore consider a spatially flat 
universe filled by dust matter, radiation, scalar field and a (negative) cosmolog- 
ical constant term. The Friedmann equations thus read : 



PM + Pr +PA + ^4'^ + V{(l)) 
1 1 



2H + 3H^ = 

where we have used natural units with 8ttG — c— 1 



3Pr.-pA + ^0'-V^(0) 



(80) 
(81) 



4.2 Matching with the data 

Notwithstanding how well motivated it is, a whatever model must be able to 
reproduce what is observed. This is particularly true for the model we are 
considering because the presence of a negative cosmological constant introduces 
a positive pressure term potentially inhibiting the cosmic speed up. Moreover, 
contrasting the model against the data offers also the possibility to constrain 
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its characteristic parameters and estimate other derived interesting quantities, 
such as qo, the transition redshift zt and the age of the universe tg. Motivated 
by these considerations, we will therefore fit our model to the dataset described 
below parametrizing the model itself with the matter density ^m, the scalar 
field quantities {^tf,,wo) and the dimensionless Hubble constant h (i.e., Hq in 
units of 100 km/s/Mpc), while we will set the radiation density parameter as 
Qr = lO^**'^ as in i93j from a median of different values reported in literature. 

4.2.1 The method and the data 

In order to constrain the model parameters we will consider several observa- 
tional test: (a)the distance modulus fi — m — M, i.e. the difference between 
the apparent and absolute magnitude of an object at redshift z, (b) the gas 
mass fraction in galaxy clusters, (c)the measurement of the baryonic acoustic 
oscillation (BAO) peak in the large scale correlation fimction at 100 Mpc 
separation detected by Eisenstein et al. ^94] using a sample of 46748 luminous 
red galaxies (LRG) selected from the SDSS Main Sample [95 , (d) the shift pa- 
rameter [55] |£| and we maximize the following likelihood taking into account the 
above test: 

(82) 

where p = (Omj ^<t>, Wq, h) denotes the set of model parameters and the pseudo- 
merit function reads : 



C oc exp 



N 

E 



gas,t 



Aip) - 0.474 
0.017 



7^(p) - 1.70 
003 



h-0.72" 
0.08 , 



(83) 



4.2.2 Results 

Table 1 shows the best fit model parameters, median values and 1 and 2cr ranges 
for the parameters (ilj\/, JIa, wq, h, fl^). 

FigslT] shows how well our best fit model reproduce the data on the SNela 
Hubble diagram and gas mass fraction. The best fit model is in quite good 
agreement with both the SNela and gas data. Actually, the values are 
respectively 206 and 48 to be contrasted with the number of datapoints, being 
192 and 42 respectively. Besides the predicted values for the acoustic peak and 
shift parameters are in satisfactory agreement with the observed ones: 

A = 0.45 , n = 1.67 . (84) 

complete discusion about this observational test and the quinstant model can be found 

in [39] 
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Figure 7: (Color online) (a) Best fit curve superimposed to the data on the SNela 
Hubble diagram, (b) Best fit curve superimposed to the data on the gas mass 
fraction. Note that the theoretical curve plots indeed fgas{z)'x{h/0.7)^'^ with h 
set to its best fit value (from \39\j ). 
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Table 10: Best fit {bf) and median {med) values and la and 2a ranges of the 
parameters (17 a/, ^^a, wqi h, fl^)) as obtained from the likelihood analysis (from 
[39]). 



Par 


bf 


med 


la 


2a 




0.283 


0.307 


(0.272,0.352) 


(0.246,0.410) 




-0.072 


-0.298 


(-0.54,-0.11) 


(-0.92,-0.02) 


Wq 


-0.72 


-0.67 


(-0.74,-0.60) 


(-0.79,-0.53) 


h 


0.632 


0.620 


(0.588,0.654) 


(0.554,0.692) 




0.789 


0.989 


(0.799,1.226) 


(0.700,1.574) 



Because of these results, we can therefore conclude that including a negative 
A leads to a model still in agreement with the data so that this approach to 
halting eternal acceleration is a viable one from an observational point of view. 


Another interesting tools to study the viability of a dark energy model is the 
point of view of structure formation. This kind of analysis can break beetwen 
models with similar prediction from the cosmic expansion history, in this sense 
the growth of the large scale structure in the universe provide and important 
companion test. Following this line in 97J the authors showed that the quin- 
stant model makes reasonable predictions for the formation of linear large scale 
structure of the Universe but it fails in the non linear regime because of the 
density contrast at virialisation increase with the value of virialisation redshif. 

Concerning the predictions of the cluster abundances, the quinstant model 
is capable of reproducing the results of the other models in a satisfactory way 
backwards in time up to redshifts a bit larger than z = 1 for the three range of 
mass values Then, it shows abrupt peaks of structure formation, in a serious 
departure of the hierarchical model for large scale structure. This seems to be 
caused by the unusual equation of state of quinstant dark energy, which behaves 
as stiff matter for redshifts a bit larger than one. This would result in enhanced 
accretion of the forming structures, both because of gravitational and viscous 
forces. 

5 Exponential Quinstant: Phase space analysis 

In this section we will investigate, from the dynamical systems viewpoint the 
quinstant dark energy model with expoenential potential V{(j>) = Vbe^''^. We 
do not consider radiation fluid here but a background of a perfect fluid with 
equation of state w = 7 — 1. The cosmological equations reads: 

^see 1391 for a further discussion about the observational resuhs and implications 
*the same behavior is obtained for other mass range |97| 
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H^-l4>^-lV{<j,)-\pM~^^~^, fc = -1,0,1, (85) 

H = -H'-- + iF(0) + A _ 1 (3^ _ 2) pM, (86) 

Pm = -SjHpM, (87) 

+ 3H<j) + = 0. (88) 

Our purpose is to re-express the former equations as an autonomous dynam- 
ical system. 



5.1 Flat FRW case 

5.1.1 Normalization, state space, and dynamical system. 

In order to get a first order autonomous system of ordinary differential equations 
(ODEs) is is convenient introduce normalized variables an a new convenient 
(monotonic) time variable. If the phase space is compact the flow of the system 
admits both past and future attractors. Let us introduce the normalization 
factor D = — A/3 > 0, the new time variable dr = Ddt, and the phase 

space variables: 

" y.^.n^" (89) 



The former variables lies in the compact phase space 

{{x,y,z):x' <l, v>Q, -l<H<l]. 

The variables ([55)) satisfy the ODEs (the prime denoting derivative with 
respect to r): 

x' = -~^xn{{j-2)x^ + {y^-l)j + 2) 

^A2/^ (90) 
y' = ly(^^-H{h-2)x^ + {y'-l)^)y (91) 

n' ^ -'^{n'^l){i^~2)x' + {y^-l)^) (92) 

For convenience, let us express some cosmological quantities in terms of the 
variables ([55]). The deceleration parameter is explicitly 



x2 (2 - 7) + (1 - 2/2) ^• 



V? 



(93) 



the fractional energy density of the scalar field is 
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and the 'effective' EoS parameter is given by 

_Ptot _ hi>^ -V[<t>) + {1-1)PM-K _ (2-7)a;2 + (l~y^)7 

(95) 

5.1.2 Form invariance under coordinate transformations. 

The system (|90ll92p is form invariant under the coordinate transformation and 
time reversal 

(r, X, y, n) (-T, -X, y, -H) . (96) 

Thus, it is sufficient to discuss the behaviour in one part of the phase space, 
the dynamics in the other part being obtained via the transformation (|96p . 
Observe also, that equations (|90II92[) are form invariant under the coordinate 
transformation y — s- — y. Then, (j90ll92p is form invariant under its composition 
with (jM)) . 

From equation ([^^ follows that H — ±1 are invariant sets of the flow. From 
equation (pij) follows that the sign of y is invariant. 

5.1.3 Monotonia functions. 

Let be defined ^ 

Z{x,y,n)=(^^^^ (97) 

in the invariant set 

{ix,y,n) ■.x^ + y'^<l,y>0,-l<n<l}. 
Then, Z is monotonic decreasing in S since 

Z' = VZ-f = -6Z {x^ (2 - 7) + (1 - ^) <o 

in S. The existence of this monotonic allows to state that there can be no periodic 
orbits or recurrent orbits in the interior of the phase space. Furthermore, it is 
possible to obtain global results. The range of Z is the semi-interval (0,+oo), 
and Z ^ as H ^ —1 (since H is bounded) and Z +oo as TL 1. By 
applying the Monotonicity Principle (theorem 4.12 [5l]) we find that, for all 
p G S, the past asymptotic attractor of p (the a-limit) belongs to = 1 and 
the future asymptotic attractor of p (the cj-limit) belongs to 7i = — 1. 

5.1.4 Local analysis of critical points. 

The system (|90II91I) admits ten critical points with the labels with i = 1 ... 5. 
In table [TT] we offer some partial information about the location, conditions for 
existence and some additional properties of them. All the critical points satisfy 
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Ti = ±1. In other words, they are solutions with H — zLD (i.e. with H ±00). 
If sign(H) = — 1 the associated solutions ends in a collapse (since H < 0), 
whereas, if sign(7i) = 1 we have ever expanding cosmological solutions. The 
expected cosmological behavior of our model is that the attractor solutions rep- 
resent collapsing solutions due the negative value of the cosmological constant. 

Table 11: Location and existence conditions of the critical points of the system 



Label 


Coordinates: {x,y,Ti.) 


Existence 


Pt 


(-i,o,±i) 


AU A 


Pt 


(0,0, ±1) 


AU A 


-'3 


(1,0,±1) 


All A 


Pt 




-V6 < A < V6 


Pt 




7 = 0, A^O 
0<7<2, |A| > 



Now, let us make some comments about the cosmological solutions associated 
to these critical points. 

The critical points and represent stiff-matter solutions which are 
associated with massless scalar field cosmologies (the kinetic energy density of 
the scalar field dominated against the potential energy density). In the former 
case the scalar field is a monotonic decreasing function of t (since its time- 
derivative is negative) . In the last case the scalar field is an increasing function 
of t since its time-derivative is positive. These solutions are always decelerated. 
The critical points P^ represent a fiat FRW solution fuelled by perfect fluid. 
They represent accelerating solutions for 7 < | . The critical points P^ represent 
solutions dominated by the scalar field {Q^ = 1, and H ±00) which are 
accelerating if A'^ < 2. Our models does not devoid of scaling phases: the critical 
points pt are such that neither the scalar field nor the perfect fluid dominates 
the evolution. There flm/^(i> = const., and 7^ — 7. 

Before proceed to make some numerical experiments let us discuss some 
aspects concerning the symmetry Observe that the critical points P^ , 

P^ , P^ and P^ are related by the transformation ([M)) with Pj^, P^, P^ and 
P^ respectively. In order to analyze the local stability of P^ , P^, Pj^ , P^, P^ 
it is sufficient analyse the local stability of P^ , P2 , Pi , P4 , P,f respectively, 
and then, infer the stability of the points in the "positive" branch by using ([96]) . 



41 



Table 12: Eigenvalues, and dynamical character of the fixed points of (|90ll92p . 
We use the notation A = (2 - "f)i2Aj'^ + X'^{2- 97)). 





l_j i ^ L, i i V cli LI 




1 IT7ti ci Tn ir'Q f'tt r a i^r 








nonhyperbolic if 7 = 2 or A = — \/6; 


Pi 


-6,-3- yfA,-3(2- 


7) 


stable (node) if A > — V6 and 7 7^ 2; 
saddle, otherwise. 


P2 


1(2 -7), -37, -f 




nonhyperbolic if 7 = or 7 = 2; 
saddle, otherwise. 








nonhyperbolic if 7 = 2 or A = \/6; 


Pi 


-6,-3+ y/fA, -3(2 - 


7) 


stable (node) if A < \/6 and 7^2; 
saddle, otherwise. 


P4 


-A2,i(6-A2),-A2 + 


37 


nonhyperbolic if A = or A"^ — 37; 
saddle, otherwise. 


Pi 


-37,|(2-7±iV^ 





nonhyperbolic if 7 = or A'^ = 37; 
saddle, otherwise. 



In table [12] we offer partial information about the dynamical character of the 
critical points corresponding to the "negative" branch. 

The critical point is nonhyperbolic provided 7 = 2 or A = — \/6- it is 
a stable node (future attractor) provided A > — \/6 and 7 7^ 2. It is a saddle 
otherwise with a 2D stable manifold and a ID unstable manifold tangent to 
the y-axis. is nonhyperbolic provided 7 = or 7 = 2. It is a saddle point 
otherwise with a 2D stable manifold and a ID unstable manifold tangent to the 
X-axis. The critical point P^ is nonhyperbolic provided 7 = 2 or A = ^/6. it 
is a stable node (future attractor) provided A < \/6 and 7 7^ 2. It is a saddle 
otherwise with a 2D stable manifold and a ID unstable manifold tangent to 
the y-axis. P^ is nonhyperbolic if A^ € {0,87,6}. Saddle otherwise, with a 
2D stable manifold provided A^ > 87 or ID if A^ < 87. The critical point Pr~ 
is nonhyperbolic if 7 = or A^ = 87. It is a saddle point, otherwise, with 
a 2D unstable manifold provided < 7 < 2 and A^ > 87. P^ {Pi) is the 
past attractor in the invariant set 7Y = — 1 provided < 7 < 2, A^ < 87 
(0<7<2, A^ >37). 

In table [HI where we present a summary of attractors (both past and future) 
for the quinstant model with fc = 0. 



In figure [8] we show some orbits in the phase space for the values A = — 




and 7=1. For this choice A^ < 87 and —\/& < A < v^. Thus the critical points 
P^ do not exist. By the linear analysis (see table [T2|) we find that the critical 
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Figure 8: (Color online) Some orbits of li9(M91]) in the phase space for the values 

A — o,''^^, 7=1. The critical points (resp. ) and P^ (resp. P^ ) are 

the future (resp. past) asymptotic attractors, P^ (resp. P^ ) having a stronger 
attracting (resp. unstable) manifold tangent to the y-axis. P^ , acts as a local 
attractor in the invariant set Ti. = 1 and P^ , acts as the local source for the 
invariant set Ti = —1 (they are, however, saddle points for the 3D dynamics) . 
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Table 13: Some properties of the critical points of the system (|90ll92p 



Label 


Deceleration q 






Pt 


2 


1 


1 


Pt 


-1 + 22 





-1+7 




2 


1 


1 


Pt 


\2 

-1 + — 

^ 1^ 2 


1 


-1 + — 


Pi 


^ 2 


37 
A2 


-1+7 



Table 14: Summary of attractors for for the quinstant model with fc = (system 
(Mil)). 



Restrictions 




Past attractor 




Future attractor 




Pi 


if < 7 < 2, < 37 


^3" 


if A < Ve, 7 7^ 2 


H=-l 


P5 


< 7 < 2, A2 > 37 


Pi 


if A > -\/6, 7 7^ 2 




P^ 


if A > -\/6, 7 7^ 2 






-1< H < 1 


P^ 


if A < V6, 7 7^ 2 




as above 








Pi 


if < 7 < 2, A2 < 37 






as above 


P,^ 


< 7 < 2, A2 > 37 



points and have a 3-dimensional stable manifold, having a stronger 
attracting manifold tangent to the y-axis (see figure [8|), i.e., two global future 
attractors might coexist (bistability). The critical points P^ and P^ are local 
sources in the invariant set = 1. Numerical inspection suggest and analytical 
results confirm that they are also global sources, P^ having a stronger unsta- 
ble direction tangent to y-axis. The critical points P^ acts locally as saddles. 
For P2 (resp. ^2*") the stable (resp. unstable) manifold is 2-dimensional and 
tangent to the y-H plane. There are orbits (corresponding to exact cosmo- 
logical solutions) connecting P^2 3 with Pj~2 3 (recoUapse occurs). The critical 
point P/, with coordinates (1/2, \/3/2, 1), have eigenvalues —2.25, —1.5, 1.5 act- 



44 



ing as a local attractor in the invariant set H. ^ 1 and , with coordinates 
(— 1/2, '\/3/2, 1), and eigenvalues 2.25,1.5,-1.5 is the local source for the in- 
variant set H — —1. They are saddle points for the 3D dynamics (see figure 

5.1.5 Bifurcations. 

The critical points (^4^, -P4 ) reduce to {P^ , -P3 ) as A ^ (VG)^. The critical 
points (^4*^, P4 ) reduce to {P^ , -Pi ) as A ^ (— V6)^- The critical points P^ 
reduce to P^ as 7 ^ 0+. On the other hand, P^ reduce to P^ as A — > (v^)^ 
or A ^ (— Y^37)^. For these values of the parameters a bifurcation arises. 

5.1.6 Typical behavior 

Once the attractors have been identified one can give a quantitative description 
of the physical behavior of a typical flat quinstant cosmology. 

For example, for ever expanding cosmologies with 7J > 0, — > +00, Ti. — I, 
i.e., the standard expanding cosmology near the big-bang, a typical model be- 
haves like a massless scalar field (kinetic dominated energy density) represented 
by P^ or P^ provided A > —a/6, 7 ^ 2 or A < -s/G, 7^2, respectively. This 
types of solutions might coexist in the same phase space. The late time dy- 
namics in 7i = 1 is given by either a scalar field dominated solution [Q^ — > 1) 
represented by P/ or by a scahng solution (Slm/170 = 0{1)) represented by 
P^ provided A^ < 87 or A^ > 87, respectively. For finite values of H, i.e., 
— 1 < 7i < 1, the early time dynamics is the same as in the previous case but 
there are subtle differences with respect the late time dynamics. In fact, in the 
invariant set —1 < H < 1 the future attractors are P^ or P^ depending if 
A < y/6, 7 7^ 2 or A > — -\/6, 7 7^ 2. If |A| < VQ the system is bistable. Such 
solutions represent contracting stiff-fluid cosmologies. This means that a typi- 
cal quinstant cosmologies allows the collapse of matter when the time evolves. 
For contracting cosmologies with {H < 0, H ^ oo,H = —1), i.e., the standard 
contracting model near the initial singularity, the late time dynamics is the 
same as int the previously described case, i.e., the collapse. However, there are 
subtle differences concerning the early time dynamics. The late time dynamics 
in ?i = — 1 is given by either a scalar field dominated solution (fi^ — )■ 1) rep- 
resented by P^ or by a scaling solution (fi,„/ri0 — 0(1)) represented by P^ 
provided A^ < 87 or > 87. 

5.2 Quinstant cosmology with negative curvature 

5.2.1 Normalization, state space, and dynamical system. 

Let us consider the same normalization as in section lS.l.H i.e, the normalization 
factor D = — A/8 > and the time variable dr ~ Ddt. We will consider 

the variables (|89p augmented by the new variable z — These variables lies 
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in the compact phase space 

{{x,y, z,H) : x"^ + + z"^ < 1, y>0, z > 0, -1 < 7^ < l} . 

The variables x, y, z, and Ti. satisfy the ASODE (the prime denoting derivative 
with respect to r): 



xH (7 - 2)x' + (t/ - 1) 7 + 7 - o + 2 



3 / VGXx 
2^ I 3 



H[{^-2)x^ + {y^-l)^ + z^j-- 



n' 



\zn ((7 - 2)x^ + {y'-l)i + z' (7 - - 

1) ((7 - 2)x^ + {y'-l)j + 7? (7 - ) 



(98) 

, (99) 

(100) 
(101) 



As before, we will re-express the cosmological magnitudes of interest in terms 
of the normalized variables. 

The deceleration parameter is explicitly 



-aa/a = —1 + — 



(2-7) + (l-y^)7 + ^M|-7) 

re 



(102) 



the fractional energy density of the scalar field and curvature are given respec- 
tively by 



-\- y'^ 

and the 'effective' EoS parameter is given by 



(103) 



Plot 
Plot 



- V{dp) + (7 - l)pM - A 



1 ^2 
2 



(2-7)a:' + (l-j/')7"7^' 
(104) 



5.2.2 Form invariance under coordinate transformations. 

The system (|98m0ip is form invariant under the coordinate transformation and 
time reversal 

(r, X, y, z, H) (-r, -x, y, z, -H) . (105) 

Thus, it is sufficient to discuss the behaviour in one part of the phase space, 
the dynamics in the other part being obtained via the transformation (I105|) . 
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Observe that equations (|98m0ip are form invariant under the coordinate trans- 
formation y —y and z —z. Then, (j98m01[) is form invariant under they 
composition with (|105p . 

There are four obvious invariant sets under the flow of (|98m0ip . they are 
y = 0, z = 0, and 7i = ±1. They combination defines other invariant sets. The 
dynamics restricted to the invariant set z = is the same described in section 
15.11 It is of course of interest the analysis of the behavior of the 4D orbits near 
the invariant set z = (or the other enumerated above). We will discuss this in 
next sections. 

5.2.3 Monotonic functions. 

Let be defined 

in the invariant set 

{{x,y, z,n) : x"^ + y^ + z'^ < 1, y > 0, z > 0, -1 < 7^ < 1} . 

Then, Z is monotonic decreasing in S since 

Z' = VZ ■ f = -Z (4z2 + 6x^ (2 - 7) + 6 (1 - _ z^) 7) < 

in S. The range of Z is the semi- interval (0, +00), and Z^OasTi^— 1 (since 
TC is bounded) and Z +00 as 7Y ^ 1. By applying the Monotonicity Principle 
(theorem 4.12 fST) we find that, for all p G S, the past asymptotic attractor of 
p (the a-limit) belongs to 7i = 1 and the future asymptotic attractor of p (the 
w-limit) belongs to 7i = — 1. 

Let be defined in the same invariant set the function 

z"^ 

Z2 = ^2 = -2 (2 - 37) HZ2 (107) 

(l-x2-2/2-z2)' 

This function is monotonic in the regions H < and H > for 7 ^ §• 

The existence of monotonic functions allows to state that there can be no 
periodic orbits or recurrent orbits in the interior of the phase space. Futher- 
more, it is possible to obtain global results. From the expression Z2 we can 
immediately see that asymptotically z ^ or x'^ + y'^ + z'^ ^ 1. 

5.2.4 Local analysis of critical points. 

The system (jQBIIlOip admits fourteen critical points. We will denote the critical 
points of the system (|98m0ip located at the invariant set 2: = in the same 
way as in the table [Tl] of section 15.1.41 We submit the reader to this table 
for the conditions for their existence. In table [T7] are summarized the stability 
properties of the critical points. 

Observe that the critical points Pg^ and are related through (|105p with 
the critical points and respectively. 
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In order to analyze the local stability of Pj^ , , , , , Pg*" , it 
is sufficient analyze the local stability of Pg" , P2 , Pi , P4~ , , Pg" , P^' respec- 
tively, and then, infer the stability of the points in the "positive" branch by 
using (|105p . In table [T7] we offer a detailed analysis of the dynamical character 
of the critical points corresponding to the "negative" branch. 

The critical point P]~ is nonhyperbolic if 7 = 2 or A = — \/6- It is a local 
sink provided A > — V6, 7 7^ 0. If A < — \/6 then there exists a ID unstable 
manifold tangent to the y-axis, and a 3D stable manifold. The critical point 
P2 is nonhyperbolic if 7 = 0, 7 = | or 7 = 2. There exist always at least 
a ID unstable manifold tangent to the x axis. The unstable manifold is 2D 
provided 7 < |. In this case there exists a 2D stable manifold tangent to the 
y-TL plane. The critical point Pg" is nonhyperbolic if 7 = 2 or A = \/6. It 
is a local sink provided A < \/6, 7 7^ 0. If A > \/6 then there exists a ID 
unstable manifold tangent to the y-axis, and a 3D stable manifold. Observe 
that P]~ and Pg" coexist provided A^ < 6. In this case both are the future 
attractors of the system, they attract solutions in its basin of attraction. The 
critical point P^ ^ is hyperbolic if A € {0, ±37,±2}. If not, it is always a 
saddle point with an unstable manifold at least ID. The stable manifold is 3D 
provided < 7 < |, 2 < A^ < 6 or | < 7 < 2, 37 < A^ < 6. It is 2D 
provided < 7 < |, 37 < A^ < 2 or | < 7 < 2, 2 < A^ < 37 or ID provided 
< 7 < |, < A^ < 37 or I < 7 < 2, < A2 < 2. The critical point P5" 
is nonhyperbolic if 7 = or A^ = 37 or 7 = §• If not, it have always two 
conjugate complex eigenvalues with positive real parts (there exists at least a 
2D unstable manifold). The stable manifold is ID provided < 7 < |, A^ > 37 

or I < 7 < |, A^ < or 2D if I < 7 < 2, A2 < Thus, P~ is always 

a saddle point. 

In the negative curvature case there are two new classes of critical points: 
Mine's solutions and curvature-scaling solutions, denoted by P^ and P^ respec- 
tively. In table [T5l it is displayed the location, existence and some properties of 
them. 

Table 15: Location and existence conditions of the critical points P^ and P^ 
of the system (|98lfT0T1) . 



Label 


Coordinates: {x,y, z,Ti.) 


Existence 


P± 


(0,0,1,±1) 


AU A 


P± 




|<7<2 
and |A| > %/2 



In table fTSl where we present a summary of attractors (both past and future) 
for the quinstant model with k = —1. 



48 



Table 16: Some properties of of the critical points Pg and of the system 

dSSHlol]). 



Label 


Deceleration q 






P± 





1 





P± 






1 

3 



In figure [9] are displayed typical orbits of the system (|98m01[) in the invariant 

set 7i = 1 for the values A = ^iid 7=1. The critical points and P^ 

are the past asymptotic attractors, P^ having a stronger unstable manifold 
tangent to the y-axis. P^ , acts as a local attractor in the invariant set H = 1. 
The Milne's universe {Pq) is the local future attractor in the invariant set y = 0- 
In figure [TO] are drawn some orbits in the invariant set Ti = —00 for the same 
choice. The critical points P^ and Pg" are the future asymptotic attractors, 
Pg" having a stronger attracting manifold tangent to the y-axis. P^ acts as the 
local source for the invariant set 2: = (they are, however, saddle points for 
the 4D dynamics). The Milne's universe (Pg~) is the local past attractor in the 
invariant set y — Q. This is the better we can do numerically since the phase 
space is actually 4D. 

5.2.5 Bifurcations. 

The critical points (^P^ , P^^) reduce to {P^ , Pg") as A — > The critical 

points (P/, P4) reduce to (Pg*", Pi~) as A ^ (— -s/B)"*". The critical points P^ 
reduce to P^ as 7 ^ 0+. On the other hand, P^ reduce to P^ as A ^ iV^)^ 
or A ^ (— ^737)^. The critical points P^, P^ and P^ coincide as 7 (|)+ 
and A {V^)^ simultaneously. For these values of the parameters a bifurcation 
arises. 

5.2.6 Typical behavior. 

Once the attractors have been identified one can give a quantitative description 
of the physical behavior of a typical negatively curved quinstant cosmology. 

For example, for ever expanding cosmologies with _ff > 0, — > +00, Ti. — I, 
i.e., the standard expanding cosmology near the big-bang, a typical model be- 
haves like a massless scalar field (kinetic dominated energy density) represented 
by Pg^ or P^ provided A > — \/6, 7 7^ 2 or A < a/6, 7 7^ 2, respectively. This 
types of solutions might coexist in the same phase space. The late time dy- 
namics in 7i = 1 is given by either a scalar field dominated solution {Q^ — > 1) 
represented by P/ provided < 7 < 2, A^ < 87 or | < 7 < 2, A^ < 2; or by a 
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Table 17: Eigenvalues, and dynamical character of the fixed points of the system 
(|98m01|) . We use the notation A (2 - j){2Aj'^ + X'^(2- 9j)). 



Label 




Eigenvalues 




Dynamical character 










nonhyperbolic if 7 = 2 or A = — \/6; 


Pi 


-6,- 


3 -^1 A, -3(2-7), - 


-2 


stable (node) if A > — \/6 and 7 7^ 2; 
saddle, otherwise. 


P2 


1(2 


-7), -37,-^,1-^ 




nonhyperbolic if 7 = or 7 = 2 or 7 = | 
saddle, otherwise. 










nonhyperbolic if 7 = 2 or A = \/6; 


Ps 


-6,- 


3+y§A,-3(2-7),- 


2 


stable (node) if A < V6 and 7 7^ 2; 
saddle, otherwise. 


Pi 




(6-A2),-A2 + 37,l- 


2 


nonhyperbolic if A = or A"' =37; 
saddle, otherwise. 


P5 


-37,1 


:(2-7±i^/A),l- 


37 
2 


nonhyperbolic if 7 = or A"^ = 37 or 7 = | ; 
saddle, otherwise. 


Pe 




-2,2,-1,-2 + 37 




nonhyperbolic if 7 = |; 
saddle, otherwise. 


Pi' 


-2, 


l±y^ -3,-2 + 37 




nonhyperbolic if A'^ = 2 or 7 = | ; 
unstable, otherwise. 



scaling solution (Jlm/fi^ = 0(1)) represented by for < 7 < |,A^ > 2; or 
by a curvature scaling solution represented by provided | < 7 < 2, A^ > 2. 
For finite values of H , i.e., — 1 < 7i < 1, the early time dynamics is the same 
as in the previous case but there are subtle differences with respect the late 
time dynamics. In fact, in the invariant set —1 <TL < 1 the future attractors 
are P^ or Pf depending if A < %/6, 7 7^ 2 or A > -%/6, 7 7^ 2. If |A| < ^/6 
the system is bistable. Such solutions represent contracting stiff-fluid cosmolo- 
gies. This means that a typical quinstant negatively curved cosmologies al- 
lows the collapse of matter when the time evolves. For contracting cosmologies 
with [H < 0,H ^ 00, H = —1), i.e., the standard contracting model near 
the initial singularity, the late time dynamics is the same as int the previously 
described case, i.e., the collapse. However, there are subtle differences con- 
cerning the early time dynamics. The late time dynamics in 7i = — 1 is given 
by either a scalar field dominated solution {il^ 1) represented by P^ pro- 
vided < 7 < 2, A^ < 87 or I < 7 < 2, A^ < 2; or by a scaling solution 
(f7m/f20 — 0(1)) represented by P^ for < 7 < |, A^ > 2; or by a curvature 
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Table 18: Summary of attractors for for the quinstant model with fc = — 1 



(system (pllTOT|)). 



Restrictions 




Past attractor 




Future attractor 




Pi 


if < 7 < 2, < 37or 






n = -i 


P5 
Pi 


if 1 < 7 < 2, A2 < 2 
< 7 < |, > 87 
if 1 < 7 < 2, A^ > 2 


Pi 


" if A < V6, 7 7^ 2 
" if A > -\/6, 7 7^ 2 


-1<H<1 


P,' 
P'l 


- if A > -\/6, 7 7^ 2 

- if A < V6, 7 7^ 2 




as above 








Pt 


if < 7 < 2, A^ < 37or 






as above 


P^ 
Pr^ 


if 1 < 7 < 2, A2 < 2 
< 7 < |, A2 > 37 
if 1 < 7 < 2, A^ > 2 



scaling solution represented by Pj provided |<7<2,A^>2. 

5.3 Quinstant cosmology with positive curvature 
5.3.1 Normalization, state space, and dynamical system. 



Let us consider the normalization factor D = \J ^ A/3+ ^ > and the 
time variable dr = Ddt, and the phase space variables: 



i^4^,y^^,H^^J^X. (108) 
V6D VSD D aD ^ ' 

The former variables lies in the compact phase space 

\^{x,y,z) : a;^ < 1, y > 0, -1 < if < 1,0 < z < l| . 

The variables x, ij, 2, and H satisfy the ASODE (the prime denoting derivative 
with respect to f): 
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x' = \xH{{2-^){l-x^)-^f)-^\v\ (109) 

y' = ^-yH{{2~^)x^+^{l-f)) + ^Xxy^ (110) 

z' = i£i7(3(2-7)x2 + 3^(i_y2) _2) (111) 

H' - -\(H^-l){{^-2)x^ + {f-l)^)+z^ (112) 

As before, we will re-express the cosmological magnitudes of interest in terms 
of the normaUzed variables. 

The deceleration parameter is explicitly 



3 

q = —aa/a? — —1 + - 



f.2 



X 



(2-7) + (l-y2)^ ' 



(113) 



the fractional energy density of the scalar field and curvature are given respec- 
tively by 

=^l±i!;r!^ = 4l (114) 
and the 'effective' EoS parameter is given by 

_ Ptot _ W - ^('^) + (7 - 1)PM - A _ ^ ^ (2 - 7) + (1 - y") 7 

(115) 

5.3.2 Form invariance under coordinate transformations. 

The system (|109mi2"]) is form invariant under the coordinate transformation and 
time reversal 

{t, X, y, z, (^-T, -X, y, z, -H^j . (116) 

Thus, it is sufhcient to discuss the behavior in one part of the phase space, the 
dynamics in the other part being obtained via the transformation pi6p . Observe 
that equations (|109mi2")l are form invariant under the coordinate transformation 
y —y and z ~z. Then, p09mi2")l is form invariant under they composition 
with (frT6| . 

There are two obvious invariant sets under the flow of (|109mi2p . they are 
y = 0, z = 0. 

From equation (|112p we can immediately see that the surfaces H = ±1 are 
not invariant provided z ^ 0. In fact, the surfaces H ~ ±1 act as membranes 
(that can be crossed). This follows from the fact that H'\h=±i — z^ > Q for 
z 7^ 0. Observe that if initially z > 0, then, from equation (|llip . follows that 
the sign of z is invariant. Only if z = 0, the surfaces H — ±1 could be invariant. 
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5.3.3 Monotonic functions 



Let be defined in the invariant set 

^{x,y,z,H) : x"^ < I, y > 0, < z < 1, -1 < 7? < l| , 
the function 

Z=- f^^, Z' = -2(2-37)H^ (117) 

[1 — — y^ ) 

This function is monotonic in the regions H < and if > for 7 7^ |. 

The existence of this monotonic function allows to state that there can be 
no periodic orbits or recurrent orbits in the interior of the phase space. Fur- 
thermore, it is possible to obtain global results. From the expression Z we can 
immediately see that asymptotically z ^ or x^ + y^ ^ I. 

5.3.4 Local analysis of critical points. 



Table 19: Location and existence conditions of the critical points of the system 

(nnmn) 



Label 


Coordinates: (x,y,z,H) 


Existence 


Pt 


(-1,0,0,±1) 


All A 




(0,0,0,±1) 


All A 


^3 


(1,0,0,±1) 


All A 


Pt 




-\/6 < A < \/6 


Pi 


(tv^?,V^\/^,o,±i) 


7 = 0, A ^ 
0<7<2, |A|> 








Pe 


(x\Q, ^1 (7+(2-7)a;*2),o) 


0<7<|, l-1<v/^ 


P± 
^7 




-%/2 < A < %/2 
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Table 20: Eigenvalues, and dynamical character of the fixed points of the 
system (|98m0ip . We use the notation A = (2 - j)i2Aj'^ + X'^(2 - 9j)) and 

M=\/i((2^*'-l)7(2" 37) + 24a;*2). 



Label 


Eigenvalues 




Dynamical character 










nonhyperbolic if 7 = 2 or A = — %/6; 




Pi 


-6,-3- y/fA, -3(2-7), - 


2 


stable (node) if A > — \/6 and 7 7^ 2; 
saddle, otherwise. 




Pi 


|(2-7),-37,-f 




nonhyperbolic if 7 = or 7 = 2 or 7 = 
saddle otherwise. 


2 

3 








nonhyperbolic if 7 = 2 or A = \/6; 




h- 


-6,-3+y|A,-3(2-7),- 


2 


stable (node) if A < \/6 and 7 7^ 2; 
saddle, otherwise. 




Pi 


-A2,i(6-A2),-A2 + 37,1- 


2 


nonhyperbolic if A = or A'^ 37 or A'^ = 
saddle otherwise. 




P5- 


-37,f (2-7±iVA),l- 


37 
2 


nonhyperbolic if 7 = or A'^ = 87 or 7 = 
saddle, otherwise. 


2. 

3' 


Pe 






nonhyperbolic 




P7- 


. Ai>_ A±\/8-3A^ (37-2)A 




nonhyperbolic if A = or A'^ = 2; or 7 = 
unstable (saddle), otherwise. 


2. 

3' 



Using pl6l) , it is possible to infer the local stability of , , P3+ , P4+ , P^ , P7+ 
from the local stability of ^3", P^^ , P{~ , P^ ,P^ , Pj . Thus, we will analyze only 
the critical points in the "negative" branch. 

The critical points P^ to ig have similar properties, dynamical character 
as unhatted ones in section 15.21 and the same physical interpretation of the 
unhatted ones characterized in section 15.11 Thus, we will not comment about 
its stability in detail. For completeness in table [^D] are summarize the existence 
conditions and stability properties. We will submit the interested reader to 
previous sections for more details. 

For the positive curvature model there are new critical points: the curve of 
nonhyperbolic critical points represent the Einstein's static universe. Spe- 
cial critical points of this family are those with the choice = (it exists 

provided < 7 < |) and x* — zfc-y/ '^'-^"^ ^ (it exists provided < 7 < |), and 

the curvature-scaling solution P^ . Concerning the stability of curvature-scaling 
solutions we have that is always a saddle point (at least one eigenvalues has 
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negative real part, and the others are of different sign). Its stable manifold is 
ID provided -\/2 < A < 0, < 7|; or 2D provided 0<A<V2, |<7<2or 
-v^ < A < 0, I < 7 < 2; or 3D if < A < V2, < 7 < |. 

5.3.5 Bifurcations. 

The critical points ^-P/, ^4^) reduce to (^P^ i -^3 ) as A ^ (a/6)~- The critical 
points (^^4+, P4 ^ reduce to (^P+ , ^ as A ^ (-\/6)+. The critical points 
reduce to as 7 ^ 0+. On the other hand, P^ reduce to P^ as A ^ iV^)^ 
or A ^ (-V37)"- {P^,Pi) reduce to (P/: A") as A ^ (V2)~ and {P^,P^) 
reduce to {P^,P^) as A ^ (— a/2)^. For these values of the parameters a 
bifurcation arises. 

5.3.6 Typical behavior. 

As a consequence that H = ±1 are not invariant sets, the determination of past 
and future attractors is more simpler. If A < — V6, 7 7^ 2, then P^ is the past 
attractor and Pg" is the future attractor. If — \/6 < A < -\/6, 7 7^ 2 the past 
attractors are both / g and P^ and the future attractors are both P^ and Pg" . 
Finally, if A > \/6, 7 7^ 2 then P^ is the past attractor and P^ is the future 
attractor. In any case the Universe evolves from a stiff regime to a stiff regime 
by crossing the value H = 0, allowing the collapse of the Universe. 

In table l^ where we present a summary of attractors (both past and future) 
for the quinstant model with k = 1. 

Table 21: Summary of attractors for for the quinstant model with fc = 1 (system 
(11091111211 'I for z > 0. 



Past attractor 


Future attractor 


P^ if A > -\/6, 7 7^ 2 


P3" if A < \/6, 7 7^ 2 


P+ if A < \/6, 7 7^ 2 


Pf if A > -\/6, 7 7^ 2 



6 Observational Test and Dynamical Systems: 
The Interplay 

Dynamical systems techniques by one way and Observational test by the other 
are strongly enough to discriminate among the wide variety of dark energy 
models nowadays under investigation. The first one is more mathematical in 
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character: they can be used to select the better behaved models, with appro- 
priate attractors in the past and future. These tools to analyse and interpreted 
results has gained a lot of attention in recent years. The techniques are so 
powerful when we want to investigate the asymptotic (and even, intermediate) 
behavior of models. The other class of tools is rather physical: they can be used 
of astrophysical observations to crack the degeneracy of classes of dark energy 
models. In the interplay, both serve to constraint the free parameters of the 
models under consideration. In our case, quintom and quinstant dark energy, 
with flat and curved geometry. 

Quintom Dark Energy Paradigm 

The model of quintom, which is mainly favored by current SNIa only, needs 
to be confronted with other observations in the framework of concordance cos- 
mology. Since SNIa offer the only direct detection of DE, this model is the 
most promising to be distinguished from the cosmological constant and other 
dynamical DE models which do not get across —1, by future SNIa projects on 
the low redshift (for illustrations see e.g. [44]). 

From the dynamical systems viewpoint we have obtained further results in 
support of the previous results in [TH [131 [TB]. For negative curvature mod- 
els, we have devised two dynamical systems adapted to the study of expanding 
(e = signff > 0) and contracting (e — signiJ < 0) models. Also, we have de- 
vised another dynamical system well suited for investigating positive curvature 
models. We have characterized the critical points of each system and interpreted 
the cosmological solutions associated. By devising well defined monotonic func- 
tions we were able to get global results for ever expanding and contracting 
models (for both negative and positive curvature models). We have reviewed 
the results concerning the flat case. It is known that, for flat ever expanding 
models the attractor will be the matter scaling solution [16] . If matter scaling 
solutions do not exist, the attractor will be phantom {w < —1) or de Sitter 
{w = —1) like. This is a difference with respect to the results in [T^ and [T3] . 
It was proved there, that the attractor solutions are de Sitter-like, unless some 
trajectories cross, transiently, the w = — 1 boundary to become even smaller 
before ending in a de Sitter phase. 

The new results we survey here are as follows: 

For negative-curvature ever-expanding models (e — sign H > 0) we have ob- 
tained the existence of scaling curvature attractors (without matter) (provided 
5 < i). The attractor solution will be dominated by DE whenever its existence 
precludes the existence of scaling curvature attractors. These solutions can be: 
phantom-like {w < —1), de Sitter- like {w = —1), or quintessence-like. This is 
a difference with respect the situation in [16j . We must notice, however, that 
if we consider other values for 7, other than 7 = 1, then the attractor of the 
system can be the matter scaling solution. This is the case ifO<7< |, 6 > ^. 
Under the above conditions on the parameters DM mimics DE. For contracting 
models (e — signiJ < 0), the attractor will be a MSF solution that mimics a 
stiff fluid. Towards the past, the typical situation is the reverse of the former 
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described. For positive curvature (closed) models, we have obtained conditions 
under which there is an orbit of type -K . This represents a cosmological 

solution starting in a ending towards a singularity describe a MSF cosmology. 
We have obtained, also, a flat FRW solution starting in a big-bang in +F y 
recolapsing in a "big-crunch" in _F. We have illustrated this results by means 
of numerical integrations of the ASODE describing this cosmological model. We 
have obtained conditions for the existence of global attractors. We have offered, 
here, only a simplified qualitative analysis (as a difference of the mathematical 
analysis in [SB] pages 69-73). However, our study have relevance by its own right, 
and can be considered in some way as a complement of the former since we have 
added a phantom field in the dynamics. We must to restate, however, that our 
analysis is not as detailed as in that reference. But its is suffice to illustrated our 
goals. The qualitative analysis in multi- scalar field (coventional) cosmologies 
with exponential potentials (in the context of assisted inflation) was done in the 
same reference, section VII, and in |89l|90], particularly for two fields. They do 
not consider phantom field as we do here. Our monotonic functions were able 
to discard the existence of periodic orbits, homoclinic orbits, or recurrent orbits. 

Quinstant Dark Energy Paradigm 



From the stability analysis of all studied models of (exponential) quinstant 
dark energy ([35], [SH]), the typical behavior, irrespective the curvature choice, is 
the evolution from an stiff-regime near the past, to a stiff-regime in the far future. 
Besides, from our dynamical analysis it seems that for negatively curvature and 
flat models, the model shows the divergence of the Hubble parameter {H) in 
the asymptotic regimes. There are other cosmological models of composite dark 
energy having stiff-matter domination as an attractor in the past, but usually 
they would not be global attractors, but local. This is the case of several models 
of quintom dark energy. From the structure formation we see that QDE makes 
reasonable predictions for the formation of linear large scale structure of the 
Universe. It reproduces reasonably well the non-linear structures from today 
up to redshifts a bit larger than one, but fails to reproduce the perturbations 
in the non-linear regime for redshifts a bit larger than one. This models are 
dynamically equivalent models of /(i?) modified gravity. It would be interesting 
to study how these f{R) models behave concerning structure formation, and 
then we would have a better understanding on how these observations would 
crack the degeneracy dark energy-/ (i?) modified gravity. 

In summary, the new results concerning quinstant dark energy are as follows: 
For the standard flat expanding cosmology near the big-bang, a typical model 
behaves like a massless scalar field (kinetic dominated energy density) and the 
late time dynamics is given by either a scalar field dominated solution (fi^ 1) 
or by a scaling solution {^jn/^(t> = 0(1)) represented by provided < 87. 
This is the standard behavior for quintessence models (without A). For finite 
values of H, the early time dynamics is the same as in the previous case but 
there are subtle differences with respect the late time dynamics. In fact, in 
this invariant set the future attractors are stiff-like (contracting) solutions with 
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H ±00. This means that a typical quinstant cosmologies allows the collapse 
of matter when the time evolves. For the standard contracting model near the 
initial singularity, the late time and early time dynamics is the reverse of the 
previously described. 

The behavior of a typical negatively curved quinstant model is similar to 
the flat situation, but not the same. The differences is in that in the limit 
H ~> 00 the late time dynamics given by cither a scalar field dominated solution 
{n^ 1) provided < 7 < 2, < 87 or | < 7 < 2, < 2; or by a scaling 
solution {Clrn/^<j> = 0(1)) for < 7 < |,A^ > 2; or by a curvature scaling 
solution provided | < 7 < 2, A^ > 2. For finite values of H, the future attractors 
are stiff-like solutions. For positive curvature models the Universe evolves from 
a stiff regime to a stiff regime by crossing the value H = 0, allowing the collapse 
of the Universe. Thus, from the dynamical view point there arc not significant 
differences between quinstant and quintom dark energy paradigms. 

Our opinion is that any dark energy model which presents a stiff-like equation 
of state in the past, during a long period of time, will predict abrupt peaks of 
structure formation, which would be the result of enhanced accretion of the 
forming structures, both because of gravitational and viscous forces. 
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